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Kin graphisches Verfahren zur Bestimmung der Zeitkonstanten 
und der Schwingungsdauer eines linearen Systems dritter Ordnung*, 


Von H. St. Stefaniak. 


1. Einleitung. In der Regeltechnik sieht man sich manchmal vor die Aufgabe gestellt, die 
Eigenfrequenzen eines linearen Systems dritter Ordnung aus den gegebenen Daten, welche oft 
Funktionen mehrerer Parameter sind, fiir eine Vielzahl von Fallen numerisch zu berechnen. 
Die bekannten Verfahren zur Bestimmung der Wurzeln einer algebraischen Gleichung dritten 
Grades sind dann entweder zu zeitraubend oder in ihrem Lisungsgang dem Problem nicht an- 
gepaBt. Im folgenden wird eine Methode angegeben, welche beide Nachteile im Rahmen der prak- 
tisch erwiinschten Genauigkeit vermeidet; dabei ist die Darstellung, ausgehend von einem ein- 
fachen mechanischen Problem, auf die Exfordernisse der Regeltechnik zugeschnitten. 


2. Beschreibung eines mechanischen Regelkreises dritter Ordnung und Aufstellung der Diffe- 
rentialgleichung. Gegeben sei etwa eine starre Platte vom Tragheitsmoment J, welche auf der 
einen Seite um die feste Achse 0 drehbar 
gelagert ist; an der anderen Kante im 
Abstande / ist sie in H, elastisch mittels 
einer Feder von der Federkonstanten c ab- 
gestiitzt (Abb. 1). Diese Feder ist ihrer- 
seits an dem freien Ende H, eines um 0’ 
drehbaren Stiitztragers befestigt. AuBer- 
dem ist in H, noch eine geschwindigkeits- 
proportionale Dampfung von der Damp- 
fungskonstanten k angebracht. Auf die 
Platte wirken beliebige zeitlich verander- : e poe 
liche Momente M= M(t)= )J/r(t) P(t), — S 
welche sie um die Achse 0 zu drehen ver- Ne 
suchen. Die Aufgabe bestehe darin, diese 
Drehungen durch geeignetes selbsttatiges 
Verschieben des Stiitztragers, welches eine 
Verlangerung oder Zusammendriickung 
der Feder bewirkt, wieder riickgaingig zu 


machen, so daf die Platte méglichst stan- WS 
dig eine bestimmte Lage, die Sollage, ein- Drocké| — 
nimmt,. Eine der méglichen Ausfiihrungen, NN 


welche die geforderte Bedingung erfiillt, 
beniitzt einen linearen astatischen Stell- 
motor, bei welchem definitionsgemaB 1 der 
Arbeitspunkt H, sich mit einer der Ver- 
schiebung des Steuerpunktes proportio- 
nalen Geschwindigkeit bewegt. Im vor- 
liegenden Falle wird der Steuerpunkt G Abb. 1, Einfacher mechanischer Regelkreis dritter Ordnung. 
mit der Platte und der Arbeitspunkt mit 

dem einen Federende in H, verbunden. Auf diese Weise erhalt man einen sog. geschlossenen 
Regelkreis, welcher von der dritten Ordnung ist, da sich dessen Verhalten, wie gezeigt wird, 
durch eine Differentialgleichung dritter Ordnung beschreiben lab. 


WN 


\ astat Stellmotor 
— pe) A3 a joenen waren 


| 
| 
| 
| 
| 
| 
| 
L 


* Diese Veréffentlichung ist eine iiberarbeitete Fassung des Abschnittes ,, Verfahren zur Lésung einer 
Gleichung dritten Grades‘: aus dem Anhang zur Habilitationsschrift des Verfassers (Techn. Hochsch. Min- 
chen 1944). ) , : 

Anmerkung bei der Korrektur: Ein einfaches Nomogramm zur Auflésung einer Gleichung dritten 
Grades findet sich auch bei R. Grammel, Ing. Arch. 12 (1941) 5. 175. 5 , 

1 R. C. Oldenbourg und H. Sartorius, Dynamik selbsttatiger Regelungen, S. 93. Miinchen u. Berlin 1944. 
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Zur Ableitung dieser Differentialgleichung betrachtet man das Momentengleichgewicht um | 


0. Mit cl2?— C, und kl? = C, muB sein 

Jg =—C,9+ Cy— Got M(t). (1) 
AuBerdem ist voraussetzungsgem4B bei der gegebenen Anordnung yp groportional -@. Das nega- 
tive Zeichen ergibt sich aus der Forderung, daf diese Anordnung von selbst eine bestimmte Lage 
der Platte, welche hier mit der Waagerechten iibereinstimmen mége, aufrecht erhalten soll. Wid 
diese um den Winkel -+ ¢ verdreht, so mu auf sie zusatzlich ein Moment in entgegengesetztem 
Sinne ausgeiibt werden. Dies geschieht durch Verschieben des Punktes H, nach oben, d. h. y 
muB bei positivem y negativ sein. Bezeichnet man den Proportionalitatsfaktor mit C,, so ergibt 
sich aus 


p= — Cyp (2) 
und (1) nach Differentiieren vnd Umordnen 
rate Ce Conn, M(t) 
amie meas conan yp eng (3) 


Die weiteren Uberlegungen kann man nach dem Vorgange von Wischnegradski\ von dem 


speziellen Ausgangssystem unabhangig machen, indem man eine Zeitkoordinatentransformation — 
einfithrt. Anschaulich ausgedriickt hei®Bt das: man mift den Bewegungsablauf statt mit einer — 


nach der Zeiteinheit 1 sec eingeteilten Uhr mit einer solchen, fiir welche die fiir das System na- 
tiirliche Zeit 


3 
af: 1 
j Rika, Pees = 4 
SS (4) 
die Zeiteinheit bildet und nennt die damit gemessene Zeit die natiirliche oderinnere Zeit t im 
Gegensatz zu der in Sekunden ausgedriickten Zeit t. Man setzt also 


t=—. (5) 


t 


Damit wird in Verbindung mit (4) aus (3) 
tt vt / M’ (t) 
Le | eR eee : 6 
as ae Si ee corre (6) 
wobei zur Abkirzung 
ss Saga ial st 8 oe | 
4=Vnez, oe B=Vay ei 


gesetzt und die Differentiation nach t mit einem Strich bezeichnet ist. Wie man aus (7) erkennt, 
sind A und B zwei dimensionslose GréBen, welche man die Hauptkonstanten des geregelten 
linearen Systems dritter Ordnung nennt. A und Berweisen sichin Verbindung mit dem ZeitmaB- 
stab t,, fiir den Regelkreis dritter Ordnung als ebenso kennzeichnende GréBen, wie es das dimen- 
sionslose Dampfungsverhaltnis D* und die Kreisfrequenz m* fiir ein gedampft schwingendes 
System zweiter Ordnung sind®?. Bekanntlich 1aBt sich eine abklingende Schwingung a €° ‘cos wit 
mit Hilfe dieser GréBen in der Form ae °° cos@* yl— D**t darstellen, woraus man durch Ver- 


gleich fiir die Kreisfrequenz w¢, die reelle Abklingkonstante 9* und das logarithmische Dekre- 
ment D* unmittelbar die Beziehungen 


*. == 2 — * 2 ee rye ey 22 D* 
We QW, yl D 3 (2° =. D > Ses yi_p= (8) 
entnimmt. 


Wenn man die durch M(t) = konst. charakterisierte Bewegung betrachtet, so 1aBt sich das 
Integral von (6) in den Formen 


p= Get + Pyeht + Pelt (9a) 
oder 
p= Det + Dy, e% cos(w, t + y) (9b) 
1 a . Shen. ; 
“SUE aes Civiling. 23 (1877), S. 95; vgl. auch E. Lehr, Schwingungstechnik. Bd. II, S. 324ff. 


* E. Lehr -‘Schwingungstechnik, Bd. II, S. 139; H. Stefaniak, Z. angew. Math. Mech. 28 (1948), S. 368 
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bzw. mit (8) sinngema8B in der Form 


p= DO, etit + Dy ,e— Pt cos (wp yi —D*t+y) (9c) 


_angeben, wobei die reellen 4,,/,,i5 bzw. 


he +A hg — Ah. ; Tema oR 
yg 28 4 BS $= 9+ iw, =o, (—D + y1—Dii) (10) 
die Wurzeln der charakteristischen Gleichung dritten Grades 
F(4, B, a) =i+4 A+ BA+1=0 (11) 
& 
Fe z + 29 
mq i 
| j 
| / 
-0,10 / 
0 
Z 
| p--4£-4™" / 
hes: il / 
e st a fe 
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ne 
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Abb. 2. Linien 2 = konst. und Linie 4 = @. 


sind. Im Gegensatz zu den entsprechenden mit einem Stern hezeichneten GréKen in (8) sind nun- 
mehr J, ~ und g reine Zahlengréfen, d. h. dimensionslos. Zwischen diesen und den Koeffizienten 
A und B von (11) bestehen bekanntlich die Beziehungen 


A=—i),—(4,+/A3) =—1,—20 =—A/,+ 2 Day, (12a) 
B= dy (42+ Az) + (Aa As) = 2 Ayo + (9 + Me) = —2 Dang, + of, (12b) 
1 = —A, (A, As) =—A (+o) =—Ao;- (12c) 


Man erkennt, daB das Verhalten eines ungestérten Regelkreises dritter Ordnung durch die An- 
gabe zweier voneinander unabhangiger dimensionsloser GréSen, der Hauptkonstanten A und B, 


16* 


i) 
N 
ne 


eindeutig charakterisiert ist. Diese bilden eine zweifache Mannigfaltigkeit von Zahlen, welche 
man nach dem Vorgange von Wischnegradski in der Ebene eines (A, B)-Koordinatensystems 
darstellen kann. Jedem Punkt dieser Ebene entsprechen nach (12) bestimmte Werte der 
Wurzeln A oder 0, @,, Wy) und D und umgekehbrt. Bei der letzten Aussage ist allerdings zu be- 
achten, da man entsprechend der zweifachen Mannigfaltigkeit der méglichen Systeme nur zwei 
von diesen GréBen beliebig vorgeben daif, um einen bestimmten Punkt (A, B) und damit ein 
bestimmtes Verhalten festzulegen; die iibrigen errechnen sich aus (10) und (12 c). Wird nun die 


Abb. 3. Linien @ = konst. 


eine davon als Parameter festgehalten und durchlauft dic andere jeweils alle méglichen Werte, 
so bekommt man in der (A, B)-Ebene eine Kurvenschar, welche die Eigenschaft besitzt, daB 
langs jeder Kurve die als Parameter gewahlte GréBe konstant ist. Nachfolgend werden dafiir 
der Reihe nach 4, Q, We, Wy und D genommen und interessierende Einzelheiten aufgezeigt. Dabei 
erweist es sich mit Riicksicht sowohl auf das Zeichnen der Kurven als auch auf das anschlieBend 


zu erlauternde graphische Verfahren als zweckmaBig, die Kurven nicht nur explizit sondern auch 
in Parameterform anzugeben. 


3a. Linien/j, = konst. Die Linien, langs welcher die (reelle) Wurzel A, = konst ist, erhalt man 


aus (12a) bis (12c) durch Elimination der iibrigen Wurzeln A,, A, in der Form 


i i i i Ingenieur-Archiy 
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1 
PA 


welche natiirlich identisch mit der Ausgangsgleichung (11) ist. Der Tangens des Neigungs- 
winkels « der durch (11a) gegebenen Geraden gegen die A-Achse ist 


dB 
=tga=—A/,. (13) 


iy 


Ly 


B= —7,A fie: (lla) 


4B 


naa cae 


70 


03 


O4 


yu 


Y 5 6 Uh 8 3 
Abb. 4. Linien w, = konst. 


Die Schar (lla) besitzt eine Hinhiillende, deren Koordinaten sowohl der Bedingung (11a) 


F*(A, B,j,) = B+1,44 ; E20 (11aa) 
als auch der Bedingung : 
aF* eee nt 

a (14) 


geniigen muB. Wie man sich leicht iiberzeugt, beriihren die Geraden (lla) die Einhillende in 
Punkten, deren Koordinaten A,, B, man aus (lla) und (14) zu 


eee, i ave (14a) 
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berechnet. Hiermit ist auch die Einhiillende in Parameterform gegeben; explizit lautet deren 
Gleichung 


) 


A?B?—4(43-+ B)+18 AB—27=0. (14b) 


3b. Linien 9 = = (A, + 42) = konst. Interessiert man sich fiir den Ort der Punkte, fiir welche: 


e=> (A, + 43) = konst ist, so hat man aus (12a) bis (12c) m, zu eliminieren; man erhalt auf 
B 
1 06 [az 9 10, ] | 
00 ‘log O% OS a 08 jy ? 
tts 102 248 i 
38 
10 ie 
Se. 
Wp=50 
Y 
fox) 
eee <a 
fod 
7H35 ] 
0. 
D6 
| 
(Op) 
e 22 A Fr1-0 
8 
o AB) 
Zul 
| ete 2B i 
4g “A 
<a 
he Ue : 
16 fy 
: 
14 
2 L 2 0-100 *+3)-0 
: Fara 
Z (48-1)2-40 28 
1 
4 | 
; $i ies 
G8 ha 
06 05 AB=1 Wp=G4 
A 
0 7 if 3 4 Oy 6 Wf 8 I 
Abb. 5. Linien w, = konst. und D = konst. 
diese Weise die Hyperbeln 
me 1 
A=—h—29, B=A20—4 (15) 
1 
in Parameterform und nach Elimination von A, explizit 
J 
el, ee ee (15a) 


mit den Asymptoten 


wi A=—2¢, B=—29(A+ 20), (16) 
die sich in den Punkten B, = 0, A, = 


—2oschneiden. Die Mittelpunkte der Hyperbeln gleiten 
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daher beim Verandern von @ auf der A-Achse entlang. Fiir das spiter noch zu erliuternde Ver- 


fahren ist es von Wichtigkeit, daB die Abszisse A nach (15) bei konstantem A, sich linear mit 0 
andert. 


Kine ausgezeichnete unter diesen Hyperbeln ergibt sich fiir 96 = 0, nimlich 
2 a ee (17) 


Der Schnittpunkt der Geraden (11a) mit dem im ersten Quadranten liegenden Ast der Hyperbel 
(17) hat, wie man sich durch nachtragliches Einsetzen iiberzeugen kann, die Koordinaten 


A, =—A,, B=. (18) 
1 
Von Interesse ist noch B 
der Ort der Punkte, fiir | 
welche 
o= = (A,+ 4s) =A, AEG A‘B (AB!) + 188-200 ; 


ist, d.i. die Verbindungs- 
linie derjenigen Schnitt- 
punkte A,, B, der Ge- 
_raden (lla) mit den Hy- 
perbeln der Schar (15), 
welche durch die vor- 
stehende Bedingung aus 
gezeichnet sind. Es er- 
geben sich dafiir aus 
(12a) bis (12) unmittel- 
bar die Koordinaten 


explizit lautet die Glei- 
chung der Verbindungs- 
linie dieser Punkte 


Bee EY da) 


3c. Linien w, = konst 
und @)—= konst. Man er- 
halt aus (12) 


es Es eee 


7 A 3 4 OF 6 7 8 


1 
ae eae oe? 
B= o* (20) 


[ @ 2 @ | Abb. 6. Diagramm zur nomographischen Auflésung einer algebraischen Gleichung 
e 


is) 


2” + we dritten Grades. 
2 
oder explizit 
A? B?— 4( A+ B?) +- 18A B— 27+ [2 w, (A? — 3 B+ 4q;)|? = 0 (20a) 
bzw. 
LS eee (21) 
(w2)® ' (m2)? = w 


Die Beziehung (21) ist mit (11) identisch, soferne man fiir — A, = 1] setzt. 
3d. Linien D = konst. Es ist mit (12) 


eee BS ee DIA (22) 
V—A; 1 
oder explizit 
(A B—1)?—4 D?[ 43+ B'— AB (4 D? + 1) + 16 D+ — 12 D?+ 3])=0. (22a) 
Aus demselben Grunde wie oben ist es von Wichtigkeit, da sowohl A als auch B bei festgehalte- 
nem A, linear in D sind. Wie man sich leicht itherzeugen kann, erhalt man aus (22a) mit D= 0 


0) S "ahd o” Ce Ve 2 e imm. ne de € stanten Ingenieur-Archiv 
e i j uU Yr Z itkons . 8 
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i ie Kinhii Diese wird, wie schon gezeigt, von jeder 
die Hyperbel (17) und mit D = 1 die Einhiillende (14b). , wie sek 
nee 7 den Punkten A, und B, tangiert (14a); auBerdem wird sie von diesen noch in 
den Punkten A,, B, geschnitten; es erfiillen naémlich (22) mit D—1 die Koordinaten 

1 —————— 
4 ee a ee (23) 

; laa ie 


Gleichung (lla), wobei nunmehr d B,/dA, + tg a ist. 


4B 


a } 70 
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Abb. 7a. Linien Pp = konst. 


Im Interesse der Geschlossenheit und Vollstandigkeit der Darstellung sind die Kurven (11a), 
(14b), (15a), (19a), (20a), (21) und (22a) in den Abb. 2 bis 5 alle wiedergegeben, obwohl einige 


von ihnen schon an verschiedenen anderen Stellen zu finden sind}, 


4. Diagramm zur nomographischen Auflésung einer algebraischen Gleichung dritten Grades. 
Man kann aus den Abb. 2 bis 5, sofern ein bestimmtes Wertepaar A, B gegeben ist, sofort oder 
nach Interpolation alle interessierenden Werte entnehmen. Fiir die Praxis ist es nun von Vorteil, 


1 F, Emde, Tafeln elementarer Funktionen, S, 40. Leipzig 1948; W. 


(1943), S. 217 sowie Grundgesetze der Regelung, S. 89. Wolfenbiittel 1947; 
Aerodyn. Labor. Techn. Hochsch., Miinchen. 


Oppelt, Arch, Elektrotechn. 37 
H. Stefaniak, Bericht 1/1942 des 
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dafS man diese auch mit Hilfe eines einzigen vereinfachten Diagrammes (Abb. 6) exrmitteln 
kann, das teilweise zwar schon von Wischnegradski angegeben, aber von ihm noch nicht 
zu diesem Zwecke benutzt wurde. Man benotigt dafiir die Grenzhyperbel (17), welche die Ge- 
biete zwischen abklingenden (D >0 oder @< 0) und anschwellenden (D <0 oder o> 0) 
Schwingungen trennt und die Kurve (14b), auf welcher einerseits die Punkte mit einer Doppel- 
wurzel liegen, und welche andererseits die Bereiche D <1 und D> 1 scheidet. Das Verfahren 


B 
17 = t 
10 
9 He | 
<5 
8 
-4 
7 t | 
p=-3 
~B 
6 = 
20 0 
5_3B 2 
2 
5 + ee 
Be 
agit 
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3 ab | 
2 
: ik 
— SS | ee steal 
0 7 2 3 Z 5 6 , 8 g 


Abb. 7b. Linien p = konst. 


besteht darin, daB man den gegebenen Punkt P(A, B) in das Koordinatensystem eintragt, von 
diesem aus die Tangente (11a) an (14b) legt und deren Schnittpunkte P, (A,, B,) mit (14b) und 


P, (A;, B,) mit (17) bestimmt. Nach (18) ist die reelle Wurzel 1, = — A, (bzw. die Zeitkon- 
stante I, — — 1/A,= B,) oder nach (13) A, = —tga. Weiterhin entnimmt man aus (12a) bis 
(12c) in Verbindung mit (18) und aus (22), (23) und aus Abb. 6, daB 

—20=A—A,, (24a) 

(2D)!= (A — 4,) (B—B,), (24b) 
bzw. a 

pa 4-4 B=8B,_ PP (G0) 
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und 
ist. 


ausreichen. Sollte das nicht der Fall sein und aus irgendeinem Grunde der Wunsch nach einer 
weitergehenden Genauigkeit bestehen, so geht man zweckmaBig so vor, daB man einen Wurzel- 


E 425510 7 a 3 4 


\Y, V V 
qo | |r Bia 448 l WL fel) 4] 
11 /) 
10 t 
9 =p 
oD 
y 
SS 
etl. | 
8 = a 
N + 
se mee 
Sia iS 
A Sia 
Q 
oP 
rs) 


=J -4f -5f -6/ -7/ -8/ -9 


Abb. 8. Linien q = konst. 


wert verbessert und damit die iibrigen aus (12a bis 12c) errechnet. Es soll nun gezeigt werden, 
wie man mit Hilfe des Diagrammes (Abb. 6) in einfacher Weise einen Schritt naher zum Ziele 
kommt. 


0 F 1 oF 

5. Linien = = p = konst und Fi = 

Verfahren einen ersten (reellen) Naherungswert A{ ermittelt, welcher gegeniiber dem wahren Wert 

4, um 6 zu groB ausgefallen sein soll, dann wird natiirlich (11) nicht befriedigt, sondern es bleibt 
ein Rest OF, welcher sich aus Af =2,+ 6A) zu 


dF = pd’ + qds2+ 643 (25) 


= q= konst. Hat man nach dem oben geschilderten 


wo? = B, (24d) | 


Die mit Hilfe dieses Verfahrens erzielbare Genauigkeit diirfte fiir die meisten Falle der Praxis — 
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errechnet, wobei man die abkiirzenden Bezeichnungen 
p=3Ai+2 Ad,+ Bund q=3/4,4+ A (26) 
2] *s . a 
eingefiihrt hat (p bzw. q sind natiirlich identisch mit bzw. aon) . Es lassen sich genau wie 
o. . oe. 1 é 
fiir die frither berechneten Werte 0, /,, @, usw. auch hier in dem (A, B)-Koordinatensystem Kur- 
ven angeben, lings welcher entweder p oder q konstante Werte annehmen. Aus (26) und (11a) 
berechnen sich dicjenigen fiir p — konst in Parameterform zu 


dl p 2 
| 2 
Ass FF 24,-4 Aa? gee Sg Pp (27) 
oder explizit zu 
A? B? — 4(A8+ B%) + 18 A B— 27 — p?(A? 3 B— p)=0. (27a) 


Fiir q = konst ergibt sich auf dieselbe Weise mit A, als Parameter 
1 
1 


Auch hier ist bemerkenswert, da A und Blinear in p und auch in q sind, beide bei festgehaltenem 
4,, Durch Vergleich stellt man fest, daB die Koordinaten des Beritthrungspunktes (A, B,) (14a) 
mit denen von (27) identisch sind, falls man hier speziell p= 0 setzt. Das gleiche Zusammen- 
fallen zweier verschieden charakterisierter Kurven findet bei (19) und (28) fiir g= 0 statt. Man 
kann daher statt (27) und (28) unter Benutzung der Beziehungen (14a) und (19) auch schreiben 


p=(#—7)—B= 2, B=2,(4—4), (29a) 
q=A+34=A—4,= =. (29b) 


Es ist deshalb nicht notwendig, die Kurven der Scharen (27) und (28) zu zeichnen, wie esin Abb. 7a 
und 7b und Abb. 8 geschehen ist, um daraus unmittelbar oder nach Interpolation p und q zu 
entnehmen, sondern es geniigt, wenn man die durch p = 0 und q = 0 gekennzeichneten auswahlt 
und die Koordinaten des Berithrungspunktes (A, B,) bzw. des Schnittpunktes (A,, B,) feststellt ; 
gemaB (29a) und (29b) entnimmt man aus dem durch die Kurve (19a) vervollstandigten Diagramm 
(Abb. 6) die gesuchten Werte p und q. Mit der Kenntnis dieser Werte ist es nach (25) dann még- 
lich, die Korrekturen 62’ zu berechnen. Als erste Naherung erhalt man dafiir 


OF 


OMe “pe (30a) 
als zweite Naherung 
eee ete Pp 2 OF OF q 
dn = — Fe VF) feo (3 yok): (30b) 


welche geniigend genau sind, solange p und q nicht von derselben Gréfenordnung wie die zu er- 
rechnenden Korrekturen 6A; sind. In diesem Falle miiBte man von der exakten Gleichung (25) 
ausgehen, welche aber die Kenntnis des erst zu berechnenden wahren Wertes A, zur Voraussetzung 
hatte. 

Mit Hilfe des verbesserten Wertes 2{’ =, —6A, kann man natiirlich in Fortsetzung des 
Verfahrens beliebig nahe an den wahren Wert A, herankommen. Es erweist sich dabei als zweck- 
miafig, in (30a) das aus (29a) mit Hilfe des genaueren Wertes 4,’ sich ergebende p einzusetzen. 

Hat man A, geniigend genau bestimmt, so ermittelt man die tibrigen GréBen mit Hilfe von 
(12a) bis (12c) zu 


SSA _y/=1__ (A+4Y ee seep zy 

ioe Sa i: al aa eae | D ee p= 2 y Ay » ae a (31) 
ta — hE | 
Orin > A omeag SS 2 


Die auf diese Weise errechneten Wurzeln erscheinen als dimensionslose Gré®Ben. Um zu den 
dimensionsbehafteten zuriickzukehren, hat man nach (5) fiir die ZeitgréBen 


t=t,, t (32a) 
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und fiir die reziproken ZeitgréBen nach (4) sinngemaB 
Atte= 12) DORI Oe (32b) 
zu setzen. 


Damit ist das Problem der Bestimmung der charakteristischen Daten der freien Bewegung 


eines linearen Systems dritter Ordnung mit Hilfe des Diagrammes Abb. 6 in einfacher Weise 
gelost. 


(Eingegangen am 20. August 1949.) 
Anschrift des Verfassers: Dr. H. Stefaniak, Miinchen 2, Walter-v.-Dyck-Platz 1. 
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Fin einfaches Verfahren zur Berechnung des Abwindes von Tragflachen. 
Von KE. Truckenbrodt. 


1. Einleitung. Die Berechnung der induzierten Abwartsgeschwindigkeit an einem beliebigen 
Ort in der Umgebung einer tragenden Linie, ist eine Aufgabe, die durch das Biot-Savartsche Ge- 
setz grundsatzlich gelést wird. Die praktische Ausrechnung der dabei auftretenden Integrale ist 
jedoch recht mithsam. Im Schrifttum sind deshalb mehrere Verfahren angegeben worden, welche 
die Ausrechnung der Abwind-Integrale vereinfachen sollen. Zuerst veréffentlichte H, Multhopp+ 
im Anschlu8 an seine bekannte Arbeit zur Berechnung der Auftriebsverteilung? ein Verfahren 
zur Abwindermittlung, das als ein mechanisches Quadratur-Verfahren charakterisiert werden 
kann. Es liefert, nachdem die Auftriebsverteilung lings Spannweite ermitielt worden ist, durch 
weitere Quadraturen die Verteilung der Abwartsgeschwindigkeit langs Spannweite in der Wir- 
belschicht fiir einen beliebigen Abstand von der tragenden Linie. 

Spater ist dieses Multhoppsche Rechenverfahren, das sich in der praktischen Durchfiihrung 
noch als ziemlich miihsam erwies, von G. Braun und H. Scharn? durch Berechnung von Kurven- 
Tafeln weiter vereinfacht worden. 

In der vorliegenden Note soll nun ein 
Verfahren zur Berechnung des Abwindes 
mitgeteilt werden, welches das Problem 
soweit vereinfacht, dafi iiberhaupt keine 
Arbeit mehr iibrigbleibt, die iiber das hin- 
ausgeht, was bei der Berechnung der Auf- 
triebsverteilung ohnehin gerechnet werden 
mu. Allerdings ist das hier angegebene 
vereinfachte Verfahren nur anwendbar fiir 
Abstande von der tragenden Linie, die 
gréBer als die Halbspannweite sind, 
was aber fiir die meisten praktischen Falle 
ausreichend ist. 


2. Ableitung der Naherungsformel. Auf Abb. 1. Erlauterungsskizze; Zirkulationsverteilung lings Spannweite. 
Grund des Biot-Savartschen Gesetzes er- 
halt man (vgl. z. B. Multhopp a.a. O.) fiir die in einem Punkt (x, y, 2) induzierte Abwarts- 
geschwindigkeit einer geraden tragenden Linie mit der Zirkulationsverteilung J"(7j) (Abb. 1) 


mes 7), 6) sia Ty Sa ces ‘|i : | 
enpeinp s Oa te; 1]> ¢) aD =f 7 ) iG 6P4a—7 PP yVe+(E—4P+ (7 n’)2 (1) 
—1 
i‘ (€ — Cai ee _.,: g | dy’ Fr 
(G—= GIP Se Gi= 1’)? ye tL (¢ Ge (7 1)? 


Hierbei bedeuten & = x ee vei = y. 5 a Geet 2 die mit der Halbspannweite dimensionslos 
gemachten fliigelfesten Koordinaten, €,(&,7) die Lage der Wirbelflache (Abb. 2) und 
y(n) = I'(n)/Vb die dimensionslose Zirkulation. 

Wir wollen die Abwartsgeschwindigkeit in groBerer Entfernung vor bzw. hinter der tragenden 
Linie berechnen, fiir Punkte in der Nahe der Wirbelschicht. Dann ist & > (¢ —¢,)?-++ (4 — 1 ye 
Entwickelt man den Klammerausdruck in (1) in eine Reihe, dann erhalt man unter Vernachlassi- 
gung von Gliedern zweiter Ordnung 


( (¢—4) —(—7') i ae we | | 


{ (GP +G—1FP [El 2 e 
Cy 5 (1 ay) 


are 2 , fod 2) 
€—tF + —7P FP ( 
_(,_£\) GW 1 
= Aye 2 an Pr 2 (ee 
1 H. Multhopp, Luftfahrtforschung 15 (1938), S.463. ? H. Multhopp, Luftfahrtforschung 15 (1938), 
S.153. 3G. Braun und H. Scharn, Luftfahrtforschung 18 (1941), 5. 179. 
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In (1) eingesetzt, ergibt sich dann 


1 
se g 1 lh (¢ ea c,)? a (4 Fa. ME if g 1 / , 
Seka ie O=(l—p)as | Soa dy [E|8 4 [ren (3) 


Das erste, mit 1/2 multiplizierte Integral stellt den induzierten Abwindwinkel in der Fliigel- 
Querebene (€ = 0) dar: 


1 
oe! ROE AGP y 
on a I MO ty FP 4 
il 


Das zweite Integral bildet die Integration der Zirkulationsverteilung iiber der Spannweite. Nach 
Kutta- Joukowsky ist der Gesamtauf- 


b 
ere: “Wle nC) +b/2 
Fare) trieb A=oV | I(y)dy, und somit 
| OT ee iter ear SARI Feo 22 
Wirbelschicht 2=2, OC Oly 0 : b2 
: wegen A=c, =~ V? Fund mitA= — 
2 F 
in fey ma 
Ca = Af y(n’) dy’ . (5) 
Z == 
Man kann also fiir (3) schreiben 
y 
Ow (&, > ¢) 
") 
| al Op 0,7, 
Ab nl Dane | @ 
ry he Sk eee) 
x | 7 4 [é%} 2A’ 
b Zur Berechnung von a,, (0, 7, ¢) wird 
auf Abschn. 3 verwiesen. 
Aus (6) kann man im Einzelnen 
folgern: 
1. Die Abwinde weit vor dem 
Abb. 2. Das Koordinatensystem fiir die Berechnung des Abwindes. Fliigel (é > 0) sind in Spannweiten- | 


, und Héhenrichtung konstant. Sie 
hangen nur vom Gesamtauftrieb des Fliigels ab, dagegen nicht von der Auftriebsverteilung. 
2. Die Abwinde weit hinter dem Fliigel (£ < 0) sind bis auf das weniger stark ins Gewicht 
fallende konstante Glied (1/4 £2) (c,/aA) gleich dem doppelten Wert des in der Fliigel querebene 
herrschenden Abwindes. 
Fir Punkte innerhalb der Wirbelschicht (¢ = ¢,) erhalt man 


™ 


ily < G 


é a 
on (Es 4) = (1 — 6) ou (04 GH) — or GaVAL 7) 
Fiir Punkte innerhalb der Spannweite |77| < 1 stellt a,, (0, 7, ¢1) = a; (7) den nach Prandtl be- 


kannten induzierten Anstellwinkel dar. 
Die Lage der Wirbelschicht hinter dem Fliigel errechnet man nach der Formel 


g 
Cy (Ss 1”) = 


[0 (1) — ae (6's 9s SAE’ (3) 
fh (1) 
wobei «(7) den Anstellwinkel des Schnittes 7 gegen die Anstrémrichtung bedeutet. | 
. An der Hinterkante (&,) mu, damit die AbfluBbedingung erfiillt ist, «(7) — ow (En- ie. | 
sein. Es diirfte genau genug sein, wenn man das Anwachsen des Wertes «(7) — a, (E, N; C1) iiber | 
ae bis zur Stelle &’ = € annimmt. Unter Vernachlassigung des zweiten Glieder in (7) ergibt | 
ich dann | 


616.9) = — 5 [a(n) — 2a (M1 [E—Sa (Ba) 


Im Gebiet vor und seitlich vom Fligel (€ > 0, beliebig und € < 0, |y| > l)ist C, = 0 zu setzen. | 
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3. Berechnung des induzierten Abwindwinkels am Ort des Fliigels. Die wichtigste GréBe in 
den Formeln (6) und (7) stellt der induzierte Abwind am Ort des Fliigels dar, «,, (0,7, &). 
Der Abwind in der Wirbelschicht ist nach (4) 

=ier 


: 1 / (¢—C )2?— (yn — 7’)? , 
Ow» (0, 4, 6) = Lim — | y Saag ie en 
w ( 1] o1) aes | (7 ) [(é fy)? (n n'y dy (9) 


Liegt der Aufpunkt auBerhalb der Fliigelspannweite (\7,| > 1), so 1aBt sich (9) in einfacher Weise 
auswerten. Liegt jedoch der Aufpunkt innerhalb der Spannweite (\7| < 1), so erfordert die Be- 
rechnung von (9) die Ermittlung des Cauchyschen Hauptwertes. Fiir die numerische Rechnung 
verwendet man am besten die von Multhopp angegebene Quadraturformel 


Ow (0, S1) =O (179) = boy y (17r) — = bin ¥ (Hn) « (10) 


_Dabei sind 5,,, b,,, feste Zahlen, die von Multhopp ein fiir allemal berechnet wurden. Die Anderung 


des Abwindes in der Nahe der Wirbelschicht wurde schon von Multhopp untersucht und lieferte 
das einfache Ergebnis 


da,, lf =Glay 
ae (Om 61) 9 fe—2, aye): (11) 
Die Kurve des Abwindes a,,, aufgetragen iiber dem Abstand von der Witbelschicht, besitzt also 
in der Wirbelschicht eine Spitze (Abb. 5). 
Der Abwind auBerhalb der Wirbelschicht laBt sich gemaB (4) durch graphische Inte- 


gration ermitteln. Macht man die Substitution 


(Oo) G9 1 1 
= ——~— <u<—, 
C€—GP + @—nF’ Tees Ware 
so kann man (4) auch in der Form schreiben 
u(r’ =1) 
1 1 

decaf a i) 
Oe (0, 1s S) = ee In y (u) du. ( ) 


4, Nachpriifung der Naherungsformel an Hand von Beispielen. An einigen Beispielen sollen 
die Brauchbarkeit und die Grenzen der in 2 abgeleiteten Naherungsformeln untersucht werden. 

Um die exakte Integration der Gleichung (1) einfacher ausfiihren zu kénnen, kann man (1) 
auch folgendermafen schreiben: 


1 
1 7) , ie sh G & | d nd? 
= — ~~ : 1 nN. (13) 
ie) are | Emer GaP +E—hP+a—vF] | 
oe 
Vorerst werden die exakten Formeln und die Naherungsformeln abgeleitet und dann wird im 


Anschlu8 daran ein Vergleich zwischen beiden Formel!n gegeben. 


a) Konstante Zirkulationsverteilung. Nimmt man die Zirkulationsverteilung kon- 
stant iiber der Spannweite an, so ergibt sich das als ,,Hufeisenwirbel* bezeichnete Wirbelsystem. 
In Ubereinstimmung mit (5) ist 


Ca 


a 


FA a= ay tet (14) 
Die Integration von (13) liefert fiir diesen Fall 
; Pind earl 
ONE ei aoe En) 
1 (15) 
jolene Sore ln Usa g | 
Soe oo (GNSS i yV2+(€—¢,)? + (47-7? 


ni =—1 


1 Man kann nachweisen, da der Integralausdruck dem Potential des Fliigels gleichwertig ist: 
il 


ae ; 
(ind =—z f ...dn . 


=] 
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Der kiirzeren Schreibweise wegen soll darauf verzichtet werden, die Formel mit den eingesetzten | 
Grenzen hinzuschreiben. Den in der Naherungsformel (6) benétigten Abwindwinkle in der Flii-. 
gel-Querebene erhalt man, wenn man in (15) & =20tsetzt. 
Fiir den Abwind am Ort des Fliigels = 0, ¢ = ¢, findet man | 
y 


| 
SY iba tl | 
“6 () = a eT (16) 


b) Parabolische Zirkulationsverteilung. Als zweites Beispiel sei eine langs Spannweite | 
parabolische Zirkulationsverteilung gegeben. GemaB (5) ist | 


(17) 


; Cn O06 1 
ty cee BES gy SERA a eto HON 
AM as omar ered needling) 
hy iszA x Ow a Ga 
60 Ca/tA Ca/tA Ca/tA “wo Cy fH Cy fA CofA = 


50 T q 30 
exakt | exokt 
—-— Noherung —-— NMiherung 


40 


60 


30 


10 


-40 


atthe 
0 O2 OF 6 G8 40 -20 
Snee tig aug see, in der Wirbelschicht in der Abb. 4. Abwindvyerteilung in der Wirbelschicht in der 
rnung € = — 1] hinter dem Fliigel bei konstanter Entfernung & = — 1 hinter dem Fligel bei parabolischer 


Zirkulationsverteilung, 


Zirkulationsverteilung. 


Die Integration von (13) liefert in diesem Fall einen etwas umstandlichen Ausdruck, auf dessen 
Wiedergabe verzichtet werden mége. Den in der Naherungsformel (7) benétigten Abwindwinkel 
in der Fliigel- Querebene erhalt man hieraus fiir = 0 und =, zu 


iy 3) 
Bee E nny (18) 


ce) Elliptische Zirkulationsverteilung. Als drittes Beispiel sei die tiber Spannweite 
elliptische Zirkulationsverteilung angenommen. GemaB (5) ist 
' Ca on lqpaaons 
y(n’) = 5 2yl —7?. (19) 
Die geschlossene Integration fiir den Abwindwinkel fithrt auf sehr unibersichtliche Formeln. 


Im Folgenden sei lediglich die Abwindverteilung 1a a 
iden g langs der Langsachse (£,7=0,6= 
angegeben, wie sie von Glauert! mitgeteilt ist, Die exakte Formel ie (oa ap 


Paes) hae apy (20) 


4 


eee (1920). Die Grundlagen der Tragfliigel- und Luftschraubentheorie (deutsch yon H. Holl) S. 148. 
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Hierin bedeutet E(k) das vollstindige elliptische Integral zweiter Gattung mit dem Modul 
k= 1//1-+ &®. Fir die Naherungsformel gilt wegen a; = a /ax A 


Ow (60,0) = 4 [1 : : ak (20a) 


epee Aas 


d) Vergleich der Naherungsformel mit der exakten Lésung. Fiir die konstante 
Zirkulationsverteilung und die parabolische Zirkulationsverteilung sind in Abb. 3 und 4 die Werte 
ie, _ Cy 
Caf. Laz, 


LZirkulationsverreilung 


a G7 G2 G3 GY GS 


_ Abb.5. Abwindverteilung auGBerhalb der Wirbelschicht 
in der Fliigelsymmetrie-Ebene 7 = 0 und der Entfernung 
€ = —1] hinter dem Fliigel. 


\ E<0 
\ 


iN Parabische Lirkulationsverreilung 


30 


\_ iiptische Lirkulationsverteilung 


60 
Lirkulationsverteiling 


40 
— exakt 
—-—WNaherung 0 esr oe 
(é) 0 Ge oY 46 08 10 
a Abb. 7, Abwindverteilung in der Wirbelschicht in der Entfernung 
& = — 1 hinter dem Fliigel fiir einen Rechteck- und zwei Trapez- 


fliigel der Zuspitzung Z = 0,6 und Z= 0,2 und dem Seiten- 
verhiltnis A = 6. 


- des Abwindes in der Wirbelschicht im Abstand 
& = —1 nach der exakten Liésung und unserer 

“ Naherungsformel verglichen. Die Ubereinstim- 

a mung ist recht befriedigend. Nach unserer Nahe- 
Abb. 6. Abwindverteilung auf der Langsachse rung hat man also das einfache Rezept, dab 


rs 0, C= bss man die Abwindverteilung «,, (€.7)) in der Wirbel- 


schicht erhalt, indem man den Abwind «,(7) am Ort der tragenden Linie verdoppelt und den 
konstanten Betrag (1/4 £2) (c,/7 A) addiert. Der Abwind auBerhalb der Wirbelschicht ist fitr die 


Wei 
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<a 


konstante und die parabolische Auftriebsverteilung in Abb.5 dargestellt im Vergleich der exakte 
Lésung mit unserer Naherung. Auch hier ist die Ubereinstimmung recht befriedigend. | 

Abb. 6 zeigt den Abwind auf der Langsachse (7 = 0, ¢ = ¢,) fiir konstante, parabolische un¢ 
elliptische Zirkulationsverteilung. Man sieht, da® in allen Fallen fir |é| > 1 die Ubereinstim: 
mung unserer Naherung mit den exakten Lisungen recht gut ist. Aber auch fiir |é| > 0,6 liefert 
unsere Niaherung noch recht brauchbare Werte. SchlieBlich zeigt Abb. 7 noch den Vergleich de? 
exakten Rechnungen mit unserer Naherung fiir drei verschiedene Fliigelformen, namlich einem 
Rechteckfliigel und zwei Trapezfliigel der Zuspitzung 0,6 und 0,2 sémtlich vom Seitenverhaltnis 
A= 6. Indiesem Fall wurde die ,,exakte“ Zirkulationsverteilung und Abwindverteilung nach 
dem Verfahren von H Multhopp gerechnet. Auch hier ist die Ubereinstimmung befriedigend| 


5. Zusammenfassung. Die Berechnung des Auf- bzw. Abwindwinkels in einiger Entfernung 
vor bzw. hinter dem Fliigel 148t sich naiherungsweise nach der Formel (6) 


CCE C) ag (! ia) Ow (9, 1), jas FB = 


durchfiihren, die fiir |&|>1 gute Annaherung liefert. Hierin bedeutet «,,(0,7,¢) den Abwind- 
winkel in der Fliigel- Querebene (§ = 0). Innerhalb der Wirbelschicht ist er mit dem induziertem 
Anstellwinkel nach Prandtl gleichwertig und kann mit Hilfe der Multhoppschen Quadratur be- 


rechnet werden. 


(Eingegangen am 20. September 1949.) 
Anschrift des Verfassers: Dr.-Ing. E. Truckenbrodt, (20b) Braunschweig, Andreeplatz 5. 
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Ausbreitungsgesetze schwacher VerdichtungsstéBe in Gasen. 


Von R. Sauer. 


1. Einleitung. In einer fritheren an dieser Stelle veréffentlichten Arbeit! wurden nichtstatio- 
nare Gasstrémungen mit ,,schwachen‘‘ VersichtungsstéBen untersucht. Dabei wurde die Entro- 
piednderung beim VerdichtungsstoB vernachlassigt und die Stofigeschwindigkeit durch den Mit- 
telwert der Ausbreitungsgeschwindigkeiten des Schalls vor und hinter dem Stof ersetzt. In der 
vorliegenden Arbeit werden im Rahmen dieser Naherung asymptotische Gesetze fiir die Aus- 
breitung schwacher VerdichtungsstéBe in einer ungestérten Atmosphare hergeleitet, und zwar 
fiir ebene, zylindersymmetrische und kugelsymmetrische StoBwellen. Solche Gesetze wurden 
schon mehrfach, z. B. von Du Mond, Cohen, Panofsky und Deeds? und von Casson und Stanton?, 
angegeben. Neuerdings wurden sie durch Untersuchungen von Friedrichs* und Anneli Lax® ver- 
scharft. Die im folgenden entwickelte besonders einfache und unmittelbare Herleitung der asymp- 
totischen Gesetze scheint jedoch bisher lediglich im Falle der ebenen Wellen® bemerkt worden zu 
sein. 


2. Aufgabenstellung und Voraussetzungen. In einer ungestérten ruhenden Atmosphare (Druck 
Po: Dichte 09, Schallgeschwindigkeit ay) breite sich ein ebener oder zylinder- oder kugelsymmetti- 
scher Verdichtungssto8 mit nachfolgender Ver- 
diinnungswelle aus. Im Weg—Zeit-Diagramm : ™ 
(Abb. 1) wird der Strémungsbereich vorne am 
StoB von der Bahn fder StoBfront und hinten am 
Wellenfu8B von der geraden Mach-Linie m mit 
der Neigung dx/dt = a) begrenzt. Ander Stof- 4, 
bahn f springen po, 09, @) unstetig auf py-+ Pp, 0 
+ 0, a)-+ a, an der Mach-Linie m schlieBen sich 
die Zustandswerte der ungestérten Atmosphare 
und der Verdiinnungswelle stetig aneinander. 
Die Druckverteilung p,+ p(x.t) in der Ver- % 
diinnungswelle ist in Abb. 1 fiir drei aufeinander- 
folgende Zeitpunkte t—0, t=—1t,, t=—t, an- 
gegeben. : 

Es wird nun die folgende Aufgabe gestellt: | 2 
Der Drucksprung p(0) und die Druckverteilung | bio) | 
in der Verdiinnungswelle | oe p(x, 0) fitr den Abb. 1. Ausbreitung eines schwachen Verdichtungsstofes 
Zeitpunkt t— ON st vorgegeben, Gesucht ist die mit nachfolgender Verdiinnung. 
StoBbahn f und die Druckverteilung p(x, t) fiir 
beliebige Zeitpunkte t > 0 und zwar insbesondere der Drucksprung p(t) an der StoBfront und 
die Breite b(t) der Verdiinnungswelle. Zur Lésung der Aufgabe setzen wir voraus: 

a) Das Medium ist ein ideales Gas mit konstanten spezifischen Warmen Cp, Cy (% = ¢p/¢y). 

b) Naherung der schwachen VerdichtungsstéBe: Der mit der Zeit t abnehmende Entropie- 
zuwachs an der StoBfront wird durch einen konstanten Mittelwert ersetzt. Dannist die Entropie 
in der Verdiinnungswelle konstant und fiir die Ausbreitungsgeschwindigkeit D (= Neigung 
dx/dt der StoBbahn f), den Drucksprung p(t) an der StoBfront und die Nachlaufgeschwindigkeit 
v(t) unmittelbar hinter der StoBfront gelten die Naherungsgleichungen’ 


D=+ [a+ (a+ @+%)}, (1) 


R. Sauer, Ing.-Arch. 14 (1943), S. 14. s 

Du Mond-Cohen- Panofsky-Deeds, J. Acoust. Soc. America 18 (1946), S. 97. 

B. Casson — J. Stanton, Bull. American Phys. Soc. 22, Nr. 2 (1947). 

K. O. Friedrichs, Comm. on applied Mathematics 1 (1948), S. 211. 

Annelie Lax, Comm. on applied Mathematics 1 (1948), S. 247. 

K. Oswatitsch, Z. angew. Math. Mech. 29 (1949), S. 130; vgl. insbes. 5. 131—132. ; 

R. Sauer, Méthodes mathématiques de la théorie des écoulements des fluides compressibles, § 28, 
Ziff.17. Paris 1950 (im Druck). 
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Pp 4x D— a 


UY 
Po #+1 a ridge )) 
dabei ist t-+ a, @ die Ausbreitungsgeschwindigkeit des Schalls hinter der StoBfront beziiglich 

der ruhenden Atmosphare. 

c) Linearisierung: Die isentropischen Zustandsanderungen in der Verdiinnungswelle werden 
wie in der Akustik linearisiert: | 

p Q 24% a 

Po eos %—la,’ (3) 
und die Wellenformen p = p(x, t)= konst» U = U(%-t)t— konst ftir den Druck, die Nachlaufgeschwin- 
digkeit usw. durch ihre Tangenten am Wellenfu® an der Mach-Linie m, also durch gerade Linien 
ersetzt (Abb. 2). 

3. Herleitung asymptotischer Ausbreitungsgesetze. Nach Voraussetzung c) sind die Wellen- 
formen p= p(X, t)r<konst usw. durch ihre Neigung an der Mach-Linie m festgelegt. Diese Nei- 
gung ergibt sich aber leicht aus den Differentialgleichungen der nichtsstationaren Gasstrémungen. 
Hantzsche und Wendt! haben sie in anderem Zusammenhang berechnet und erhielten fiir die Nach- | 
laufgeschwindigkeit v= v(x, t) — const (Abb. 2) 


“+o t+ konst SES ire ee 
oe (4+ 1)t-+ konst-jt—— («++ 1)¢ fir o=1, (4) | 
“x +1 x+1 | 


folnee tuna —— sor 


t-Int—+ konst:t > 


2 


Dabei bedeutet g = 0 die ebenen, ¢ = | die zylindersymmetrischen und o = 2 die kugelsym- 
metrischen Wellen. Rechts sind jeweils die asymptotischen Ausdriicke fiir t—> 00 angegeben. Die 
Gleichungen (4) zeigen, wie sich die Wellenform mit | 
zunehmender Zeit abflacht. 

Nach Festlegung der Verdiinnungswelle durch d x/dv 
bleibt lediglich noch die Aufgabe, die jeweilige Lage 
des VerdichtungsstoBes zu finden. Hierfiir beniitzen 
wir die geometrische Beziehung (Abb. 2) | 


t 


und die Gleichungen (2) in Voraussetzung b). Dann 
hat man sogleich 


P ab 
Po * da 49) (>) 
und 
Rover db x+l1dv 
Abb. 2. Linearisierte Verdiinnung lle. ——! — = —— — 
swelle dt ay A Flog © (6) 


Nach Einsetzen der asymptotischen Ausdriicke der Gleichungen (4) erhalt man durch eine 
elementare Integration aus (6) die Breite b(t) der Verdiinnungswelle und hierauf aus (5) sofort den 
Drucksprung p(t) an der Sto®front. Wir stellen das Ergebnis in folgender Tabelle zusammen: 


Pp (t) pl? ld 1 

Po t(Int)'”” 
b(t) ; 1/2 1/4 1/2 
ee ~t ~ ~~ 

(0) t (Int) 


Natiirlich kann man in diesen asymptotischen Beziehungen die Zeit t durch den Weg x= aot er- 
setzen, den der WellenfuB zuriickgelegt hat. ; : 

4. Diskussion. Die in der Tabelle zusammengestellten asymptotischen Gesetze gelten fiir 
hinreichend grofe t, fiir welche die Druckspriinge so klein sind, daB die Naherung der schwachen 


1 W. Hantzsche und H. Wendt, Jahrb. Dtsch. Akad. Luftfahrf h : i 
Bericht Angew. Math. Teil UI, Abschn. 1V, Gasdynamik, S. 122. ian paerrighingier ie ee 
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StéBe (Voraussetzung b) und die Linearisierung (Voruassetzung c) keine ins Gewicht fallenden 
Fehler verursachen. In solchen t-Bereichen werden die Gesetze durch Messungen gut bestitigt. 
Sie gelten auch fiir sgezahnartige Druckverlaufe (Abb. 3) mit einer auf beiden Seiten durch einen 


VerdichtungsstoB abgeschlossenen Verdiinnungswelle. 


lassen sich aus zwei an der Mach-Linie m aneinander- 
gefiigten Dreiecksdruckverlaufen, wie wir sie in dieser 
Arbeit betrachtet haben, zusammensetzen. 

Man beachte, daf die vorliegenden Untersuchungen 
trotz der als Voraussetzung c) eingefiihrten Lineari- 
sierungsvorschriften tiber den Rahmen der Akustik 
hinausgehen. In der akustischen Naherung erhielte 
man an Stelle der StoSbahn f eine zu m parallele 
Mach-Linie, also eine zeitunabhingige Wellenbreite 
b= konst, und auch fiir das Abklingen des Druck- 
sprungs p(t) andere Gesetze, z. B. im ebenen Fall 


(c= 0) einen zeitunabhangigen Drucksprung p= konst. 


Solche sagezahnartigen Druckverlaufe 


PotP 


Abb. 3. Sagezahnartiger Druckverlauf. 


SchlieBlich sei noch darauf hingewiesen, 


daB die Absorption, die im akustischen Bereich fiir die Dampfung der Schallwellen schlieBlich 


mafgebend ist, bei der vorliegenden Untersuchung naturgemaf vernachlassigt wurde. 


(Eingegangen am 12. Oktober 1949.) 
Anschrift des Verfassers: Prof. Dr. R. Sauer, Minchen 19. Romanstrafe 34. 


Berechnung eines statisch unbestimmt gestiitzten und senkrecht zu seiner 
Ebene beliebig belasteten geschlossenen Kreisringes. 


Von W. Swida. 


1. Einleitung. Die ersten Unters uchungen dieses Problems beziehen sich auf einen gleichmaBig | 


belasteten Kreisring mit gleichen Abstanden zwischen den Stiitzen. Dieser Fall wurdevon K. 


Federhofer!, F. Diisterbehn? und S. Hefler? und am eingehendsten von G. Unold* behandelt. 


Unold hat gezeigt, da®B die Beanspruchung des Ringes von der Profilart des Querschnittes ab- | 
hangig ist. Dabei ist zwischen Ringen ,,erster Profilart*‘ (Rechteck, Rechteckrohr, Kreis) und | 
,zweiter Profilart’* (I- oder [- Querschnitt) zu unterscheiden. Die Untersuchung von G. Unold | 


wurde durch F. Stiissi> erweitert. 

Eine allgemeine Lisungsmethode fiir einen Kreisring der ersten Profilart ist fiir beliebig ver- 
teilte Belastung von C. B. Biezeno und J. J. Koch’ auf dem Wege der Integration der Differential- 
gleichung der elastischen Linie angegeben worden. Dabei wurden als Unbekannte die vertikalen 
Stiitzkrifte gewahlt und die Belastungsfunktion in Fourierreihen entwickelt. 

In der vorliegenden Arbeit wird ein Verfahren zur Berechnung eines statisch unbestimmt ge- 


stiitzten und senkrecht zu seiner Ebene beliebig belasteten Kreisringes der ersten Profilart an- | 
gegeben, ein Verfahren, das sowohl einfach, als auch geniigend allgemein ist. Dabei wird eine so- | 
zusagen statische Verwandtschaft zwischen einem Kreistrager auf mehreren Stiitzen und einem — 


geraden Durchlauftrager ausgeniitzt. Obgleich die Untersuchung des Kreistragers durch das Vor- 
handensein veranderlicher Torsionsmomente erschwert wird, laBt sie sich doch auf die Auflésung 


eines Systems dreigliedriger Gleichungen zuriickfiihren, wie es auch bei einem Durchlauf- | 
trager der Fallist. Die Auflésung des Gleichungssystems liefert die Biegemomente an den Stiitzen | 


durch die alle anderen statischen Gréfen ausgedriickt werden kénnen. 
Es wird angenommen, daB® alle 


Pais Auflager von der in Abb. | dargestell- 
wpe ten Art sind. Ein solches Auflager 
m7 @ 7 +7 


behindert die Drehung des Stiitzen- | 
querschnittes in der normalen (ra- 


dialen) Ebene, gibt ihm aber die volle 
Méglichkeit, sich in der Tangential- 
ebene NO (Abb. 2) zu drehen. Auf 
Abb. 1. Abb. 2. diese Weise hat jeder Stiitzenquer- 
schnitt nur eine Drehungsfreiheit 


und zwar in der Tangentialebene des Ringes. Der Fall eines ,,freien‘ Kreisringes, bei dem | 


jeder Stittzenquerschnitt immer zwei Drehungsfreiheiten hat, bei dem also nur der Schwer- 
punkt des Stiitzenquerschnittes unbeweglich ist (punktweise Stiitzung), ist als Gegenstand fiir 
eine gesonderte Untersuchung vorgesehen. 


2. Theoretische Grundlagen. a) Allgemeines. Fiir eine beliebige Stiitze ist 
0 =— @,. (1) 


Dabei sind @ und @, die Neigungswinkel der Tangente an die Biegelinie in der vertikalen Tan- 
gentialebene NO (Abb. 2). Diese Gleichung kann man folgendermafen darstellen: 


An Wnt ie = — Ane eee 


oder 


Nn Nn bn + bln =— (A, An). (2) 
Hier bedeuten 4,, 1), n> n> ln und uy, die Drehwinkel an der Stiitze n fiir die Bégen s, und s 


* K. Federhofer, Z. Math.u. Physik 62 (1914), S. 382 
2 F. Diisterbehn, Eisenbau 6 (1920), S. " : ig 
; - ee Beton und Eisen 29 (1930), S. 149. 
- Unold, Forschungsarbeiten auf dem Gebiete des I i i 
5 F, Stiissi, Schweiz. Bauztg. 111 (1938), $.166.. © slate 


° C. B. Biezeno und J. J. Koch, Z. angew. Math. Mech. 16 (1936 5 Si L 
RaGtzmmel, Technische Dynandk, S338, Berlin 1030. 2 He ee 


n+19 
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die als freiaufliegende Kreistrager auf zwei Stiitzen (Abb. 1) zu betrachten sind. Hierhei ist A, der 
Drehwinkel des Elementes s,,infolge der AuBenbelastung des Feldes s,; 2, ist der Drehwinkel des 
Bogens s,,, infolge seiner Aufenbelastung. Ferner bedeutet Yn baw. 7, den Drehwinkel fiir den 
Bogen s,, bzw. s, +, infolge des Biegemomentes M,, an der Stiitze n. M,ist der Drehwinkel des Ele- 
mentes s, infolge des Biegemomentes M,_; an der Stiitze n — 1 und schlieBlich ist Mn der Dreh- 
winkel des Elementes s, 4) infolge des Biegemomentes VES 

Wir ermitteln zunachst die Werte 7,, 7/n, /, und yu, fiir die Momente, sowie die Drehwinkel 
A, und A, fiir einige gebrauchliche Belastungsfiille. 

b) Die Drehwinkel jj und yw. Man betrachtet einen waagerechten Kreisbogen AB (Abb. 3). 
der auf zwei Stittzen, nach der Art der Abb. 1, frei aufgelagert und durch das Biegemoment M 
an der Stiitze B beansprucht wird. Das Moment M wirkt in der vertikalen Tangentialebene des 
rechten Stabendes. Unter der Wirkung von M entstehen an den Stiitzen die Auflagerkrifte A 
und B und die Torsionsmomente m, und m,, die in den Ebenen der Stiitzenquerschnitte wirken. 
Die Pfeile der Momentvektoren sind in der itblichen Weise nach jener Seite gerichtet gezeichnet, 
von der aus gesehen die Momente im Uhrzeigersinne drehen. 

Man ermittelt zunachst die Auflagerkrafte und die Torsionsmomente m, und m,. Die Gleich- 


gewichtsgleichungen des Ringabschnittes AB sind 

DL =A B=; 

> m, = M sina — (m, +m,) cosa —0, 

ym, =2 Arsina +(m,— m,) sina — Mcosa =0. 
Die restlichen drei Gleichgewichtsgleichungen gehen in die Identitat 
0 =0 itber. Daraus folgt 


M 
Urs rae 
m, = Ar— M-ctg2oa. 
Zur Ermittlung von A wird der Satz von der kleinsten Formanderungs- Abb. 3. 


arbeit beniitzt. Das Torsionsmoment fiir einen beliebigen Querschnitt 
(ab), der einem verdnderlichen Zentriwinkel y entspricht, ergibt sich als algebraische Summe 
der statischen Momente aller links von (ab) angreifenden Krafte in bezug auf die Tangente 


aK) zu M 
ee) M, =m, cos y — Ar (1 — cosy) = a eee 


Das Biegemoment im Querschnitt (ab) als algebraische Summe der statischen Momente derselben 


Krafte beziiglich des Radius (aQ) wird 


Msi 
M, =Arsing + m,sing = aad 


sin2a ~ 

Der positive Richtungssinn des Radius (a0) und der Tangente (aK) ist in Abb. 3 jeweils mit 
einem Pfeil bezeichnet. 
Das Biegemoment kann auch aus der Bedingung: 


bestimmt werden. 
Die potentielle Energie des Kreisbogens AB wird 
M?} ds M? ds 
V | 
ee a es 
(AB) (AB) o." 

wobei C die Torsionssteifigkeit ist. Der EinfluB der Schubkrafte kann dabei vernachlassigt 
werden. Ae = 

Auf Grund des Satzes von dem Minimum der Formanderungsarbeit ergibt sich fiir einen Ring 


mit konstantem Querschnitt : 


20 Diy 
0M 
ay I aM, 1 OMAee 1 om; Lf msg dae = 
oy hm, ids + th Mag ds zy ee ec oa 
(AB) (AB) 0 ° 
20 


- nee Ar) dg =— 5 (M—2 Ara) =0. 
C 


sin 20 C 
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Es ist also 


M 
A Pea RaM 
1 1 |. 
3) iM 
m = M(s35 za) 8) 
mie ya ete 2a), (4) | 


rs cosy 1 
M, =M(s5 cake 
Nun errechnet man den Drehwinkel 7 an der Stiitze B. Auf Grund des Satzes von Castigliano ! 
wird 


24 
ay 1 om Meany lo, Mi ae ee [MESOe ay 4 
=a b—— Des M, ay ds = P 


ako oM EJ } “sin? 20 
(AB) (AB) 0 
20a 
. 5) | 
1 cos p Te\e Mr ( Bal ( | 
ee | OE a ara Og A ie 't 
if, 
Mr a ctg 20 1 Mr Mr 
a cale ae ayes Ea (a) += fe (a). 


Um den Drehwinkel am linken Ende A zu bestimmen, laBt man an diesem Ende ein fingiertes | 


Moment M’ wirken (Abb. 4). Dann wird 
M — M’ 
27TH 


A 


Aus den Gleichgewichtsgleichungen 
(M + M’) sin « — (m, + m,) cosa =—0, 
M — M’ 
2ra 


2 rsin a +(M’— M) cosa +(m,— m,) sina —0 


ergibt sich 
M . M— M’ 
Ce oe ce 2 eee 


Man erhalt schlieBlich 
M, =m, cos y — M' sin gp — Ar (1— cos) = 


= cos p 1 off ut ; 
Were Fe m'(y, ctg 2.4 cos p— sing), 


Abb, 4. 


M, =Arsingy +m,sin y + M’ cos Sey + M’ (cos my — ctg 2a sing) , | 


sin 2a 
aM. fake | 
a OV —_ 1 is il M, ; 
- ae (ay | Me satrap | [Msn rap Preaen 
6 0 
20 20 
_ 1 { Mrsing i di cos Y 
=r sin 2a (COSP — cts 2a sing) dee | | (6) | 
i 0 | 
1 I : M 1 
— Fh) (Fy — 008 tg 2a — sin g) dp = MY gre (1 — 2avotg 2a) 
Mr/ 1 1 a ctg 2a Mr Mr 
Weg 8 = 
baa 6; i 2 sin 2a oe +e fa (a) - 


_¢) Die Drehwinkel 4 fiir gleichmaBig verteilte Belastung. Die Auflagerwider- 
stande sind A — B —gqrq (Abb. 5). Man bringt an den Schnittstellen 4 und B die fingierten 
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Momente M’ an. Aus der Gleichgewichtsbedingung 
a 
>) m, =2 M’' sina — 2mcosa +2f gr (reosy — rcosa) dy =0 
) 6 
folgt 


m = qr" (tga—a)+M’'tga. (7) 
Fiir einen beliebigen Querschnitt (Abb. 6) ist 


P 
M, =m cos y — qr? «(1 — cos 7) + fqrdé(r—rcos é) — 
U 
— M’sin » = qr’? (tg «cos p— a + — sing) + M’ (tg « cos y — sin 9), 
” 
M, = qr’ asin p +m sin p — | qr®sin Edé + M’ cos 
v 


=qr* (tg asin —1-+cos@) + M’ (tg a-sing + cos Q) 


Abb. 5. Abb. 6. Abb. 7. 


Der Drehwinkel am linken bzw. am rechten Ende wird 
SES ia ‘oes: ene ' , F 
= Flaw) wo DES! I" (te & sin p — cos y) (cosy +tg asin) rdp + 
0 


20 


= qr? (tg a: cosy — a + — sin @) (tg a cos y — sin ¢) rdp (8) 


0 


3 3 3 3 
Se —tga) +27 (atgta —3tga $3.0) = f(a) +2 f(a). 


cos? a 
Die Zahlenwerte der Funktionen f, (x), fy (a), ... fg (x) in den Formeln (5), (6) und (8) sind 
weiter unten in einer Tabelle aufgefiihrt, und zwar fiir den Winkel a gleich 45°, 30°, 22° 30’ und 
ae 
d) Die Drehwinkel / fiir eine unsymmetrische Kinzellast. Die Lage der Einzellast 
wird durch den Winkel f angegeben, wobei 6 < « ist (Abb. 7). Man Ja8t am linken Ende des 
Kreisbogens AB ein fingiertes Moment M’ wirken. Dann wird 


Py M’ 


Aus den Gleichgewichtsbedingungen 
> m, = Pr (cos (« — B) — cos a] — m, cos a — m, cosa + M’sina =0, 
ym, =2 A-rsin x — Pr [sin (a — B) + sino] +m, sina — m,sina + M’ cosa =0 


erhalt man 


Sed apres sunt 8) =F +m" : mae (9) 


sin 2a 20 


m =Pr(5 7) , M'( eee on 
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wobei 


ae sin (x — fs) ! cos (a —f) 
sing cos & 


ist. Fiir einen beliebigen Querschnitt des Abschnittes A P wid 


Mi =m, cosq — Ar(1—cosqy)— M’ sing =Pr (‘ “S Pet) Z| $M (t ctg 2a cosy — sing), 


M; =Arsing 4m, sing + M’ cosy =*" ksin y + M’ (cos gp — ctg 2a sing). 


Fiir einen beliebigen Querschnitt des Abschnittes PB wird 
Mi" =m, cosy — Ar (1— cos gy) + Pr [1 — cos (vy — f)| — M’ sin p 


it 
== Pr E cos 1 — cos ( — 8) + M’ (5— ctg 2a cosy — sing), 
My =Arsing +m,singy — Prsin(y— f) + M’ cosy = 
= Pr aae sin (py — 6) + M’ (cosg — ctg 2a sing). 


Der Drehwinkel an der linken Stiitze ergibt sich zu 


} 
F 0 
NA = oes k sin @ (cos y — ctg 2a sin gy) rdg + 
B 
ee 4) & ctg 20 cos p — sing) rd + 
2 
20 
1 — 
ay — sin (y — f)| (cos y — ctg 2.4 sin y) rdy + | 
2 | 
20 
7 k 1 : al 
are 1 ~ cos p +# — cos (p — 8) J —etg2.acos p — sing) rdp = (11) 
B 
Pr { si in (4a — k k sin 4 : 
=F hore 20 | ee B) - : oe (x A) cos] +sin p(x —§)\ 4 
Pr? { ee sin(4a—B) ka ksin4« 
4 Ee ete 2a [2 sin p A 5 | (3 Ecos A] + 


k /sin2 a ; B 3 sin (2 «4 — a 
| in? 2 | pi a—B) cos (4a—f) 
(PS —sint 20) 41 — feos p— GRD _ cep 


| 
ans 
| 
to|D 
Pa 
a. 
=} 
oD 
CaSO 


; ; Pr? TAO 
Die Beiwerte von EJ und Ge in der letzteren Gleichung sind von den Zentriwinkeln « 


und § abhangig. Fiir zwei Laststellungen und zwar fiir B —a/2 und f —2 4/3 sind diese ‘a 
werte in der Tabelle 1 angegeben. Dabei bedeutet f, (a, 8) den Beiwert von aa fiir B =F 


: es 2 
und f; (a, f) denselben Beiwert fiir f == . Ferner ist f, (x, 8) baw. f, (x, 6) AG Beiwert von 


P 
fir 6 =% baw. p=28. 


Zur Bestimmung des Drehwinkels Ag am rechten Ende wird daran ein fingiertes Moment M’ 
angebracht (Abb. 7). Dann erhalt man 


oy) cas ae ee ic 
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2 20 sin 24 20 


ae | ! mt (288 ‘). 


Py M’ sin 
my Pr. : Q 
: 2 FS SP San Shige. © 
ie k 5 8 ,{ Cos 1 
mf, = Pr cos P15, ¢o8 (p —p)| +. M ee — ) : 
II [oe : M' sin 
M;, = Pr|* sing (p— p)| 4 ar gee 
B B 
Ue | eR reesin 0 if k cos @ y cos ~ i 
oe ay eat! af 
Ap aI 2 k ainda’ ?? =| Pel 2 =) ie 26 2a rdg | 
0 0 
Ge 
I P bk. é sin ~ ve ie 
re ee ars eo 
B 


; “ 3 A cos P 1 
+4? [Ecos 9 +35 — 08 (g A & 2a zs)ré9 


—-EJ-sin 2.0.| 2 8 4 4 


272) al ka ksin 4a sin (4 « —f) of is B 
Pr | ie ; nl 7 sin B (af) cos 8] + FY. 


In der Tabelle bedeutet f; (a, 8) den Beiwert von Pr?/E Jin der letzteren Gleichung fiir 8 = a/2 
und f, (a, f) denselben Beiwert fiir 6 = 2.4/3. Ferner ist f, (a, 6) baw. f, (a, 8) der Beiwert von 
Pr?/C fir B =a/2 bzw. fiir B =2 «/3. 

Wenn die Einzellast Pin der Feldmitte wirkt, dann ist / —«, und man erhalt nach einer Reihe 
von Umformungen 


| = 
GhiGe\sin 2a 


eee TOUR ir / ote 8 1 aed as IPP 
Mde= ar gene tC aeaa Tie te) TO t+ eho). (3) 
Dabei wird das Torsionsmoment an der Tabelle 1. Funktionen f(x) und f (a, f). 
Stiitze Funktionen f(a) bzw. f(a, B) fir 
— at is 1) ; (14) a = 2/4 mo | m/8 | — #/10 
BN f(a) 0,7854 | 0,4095 | 0,2854 | 0,2211 
Die Funktionen f,(a) und f,(x) sind eben- — f2 (4) 0,1488 | 0,0320 | 0,0122 | 0,0060 
fy (x) 0.5000 | 0,2283 | 0,1517 | 0,1150 
falls in der Tabelle angegeben. es 0.1366 | 0.0285 | 0.0107 | 0.0054 
Es ist zu bemerken, da samtliche Bei- fz (ax) 0.2854 | 0,0604 | 0.0229 | 0,0112 
werte f(x) und f(a, f) fir 20 —180° ent- Ff, (a) 0.0708 | 0,0067 | 0,0014 | 0,0005 
weder unendlich gro sind, oder in eine un- i e ee Heke Rhaoe oe 
. +. G x ) o) 9 Cl ? 
bestimmte Form itbergehen. Dasselbe gilt 4G.P) 0.2255 | 0.0750 | 0.0385 | 0.0230 
fiir die Torsionsmomente an den Stiitzen. f, (a, B) 0.0530 | 0.0077 | 0,0018 | 0,0608 
In der Tabelle wird der EinfluB der Biege- —_—f, (a, f) 0.1791 | 0,0542 | 0,0269 | 0,0160 
momente durch die Beiwerte f, (a), fs (4), i ne poe Meat eee ee 
1 a, r) y) ° ’ 
-- fi (6), f(a, B) usw. mit ungetaden  (.°% | 9'0630 | 0.0090 | 0,0023 | 0,0010 
Zeiger und der EinfluB der Torsionsmomente — f(x.) 0.2260 | 0,0692 | 0,0340 | 0,0200 
durch die Beiwerte f, («), f,(),---fo(0B).fa  f,(a.B) | 0,0583 | 0,0073 | 0,0019 | 0,0008 


(x, 6) usw. mit geraden Zeiger beriicksichtigt. 
e) Dreigliedrige Elastizitatsgleichungen. Nun kann man auf Grund von (5) und (6) 
folgende allgemeine Ausdriicke fiir Drehwinkel 7,,, Yn> fn und pi, in (2) anschreiben: 


M,r Mer. 
n= Ry hi (an) +—a-S (Xn) ; 
M_r M,r 


tn = ep fi (%m41) +E Sa (Ont) » 
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Moat Mit 
Pa = ae (on) as G fa (an) ? 


ty = EEE Fa Gomi) cp a faa 
Setzt man die erhaltenen Werte von 7, Yn Mn und fu, in (2) ein und multipliziert die beiden Seite | 
dieser Gleichung mit EJ/r, so erhalt man 
M,,-1 Lf's (On) +t fa (On)] + Ma [fi (On) +A (On41) + tf (On) | 
+t fy (%n+1)] + Mn4t [fs (+1) +t fa (%n+1)] (15)| 


wobei t= EJ/Cist. Das ist die Grundgleichung zur Ermittlung der statisch unbestimmtem| 
Stiitzenbiegemomente in einem geschlossenen oder offenen Kreisringtrager. Diese dreigliedrige 
Gleichung ist ahnlich gebaut wie die Clapeyronsche Dreimomentengleichung. Die Beiwerte vom 
M,,_1, M,, und M,,,; kénnen als Vorzahlen und die rechte Seite der Gleichung als Belastungsghiedt 
bezeichnet werden. 

Wenn die Zentriwinkel gleich sind, wenn also a, =a, 11 = ist, dann vereinfacht sich dies 
Grundgleichung wesentlich: 


(M,—1 + Ma +1) [fs (a) +tfa (0)] +2 Ma Lh (0) + tf (@)] =— An +n) 


Fiir jede Stiitze eines geschlossenen Kreisringes mu eine Gleichung nach der Art (15) oder 
(16) aufgestellt werden. Die Zahl der Gleichungen muB also gleich der Anzahl] der Stiitzen sein. 
Fiir symmetrische Belastung kann natiirlich eine wesentliche Vereinfachung erzielt werden. 

Fir einen offenen Kreistrager auf n Stiitzen (Abb. 10) geniigen n — 2 Gleichungen fiir die: 
Zwischenstiitzen. 


EJ H 


| 
3. Zahlenbeispiel. Ein geschlossener Kreisring auf fiinf Stiitzen mit r= 3 m hat die in Abb. 85 
dargestellte Belastung aufzunehmen. Die Stiitzenabstande sind verschieden. Der Einfachheit 
halber wird angenommen, daft = EJ/C =1 ist. 
Wir fangen mit der Stiitze 2 an. Die Gleichung (16) fiir diese Stiitze lautet 


(M, + M3) [fs(a) + tf, (a)] +2 M,[fi (a) +f, («)] = 
EJ [Pr Pr r3 8 
=F [Eh +The) + Ez hl) +55 Jol@). 
Nach Einfithrung der Zahlenwerte erhalt man 
(M, + M,) (0,2283 + 0,0285) + 2M, (0,4095 + 0,0320) 
= — [5-3 (0,872 + 0,0099) + 2-3? (0,0604 + 0,0067)] . 
Fiir die Stiitze 3 ergibt sich nach (15) 
M, (0,2283 +-0,0285) +. M, (0,4095 + 0,0320 + .0,7854. -+ 0,1488) + M, (0,5000 + 0,1366) 
= — [2-3? (0,0604 + 0,0067) + 3,5-3 (0,2255 +. 0,0530 +.0,1791 + 0,0482)] . 
Es ist leicht, die restlichen drei Gleichungen aufzustellen: 
(M; + M;) (0,5000 + 0,1366) +2 M, (90,7854 +- 0,1488) 
= — [3,5-3 (0,2255 +.0,0530 +.0,1791 + 0,0482) + 4,5 +3 (0,2260 + 0,0583)] (Stiitze 4), | 
M, (9,5000 -+-0,1366) + M, (0,7854. + 0,4095 + 0,1488 + 0,0320) + 


+ M, (0,2283 +-0,0285) = — 4,5-3 (0,2618 + 0,0630) (Stiitze 5). 
(M; + M,) (0,2283 + 0,0285) +2 M, (0,4095 +-0,0320) =— 5-3 (0,0872 +-0,0099) (Stitze 1). 
Die Auflésung der Gleichungen liefert | 
M, =— 0,62tm, M,=—2,06tm, M, =— 2,76tm, M, =— 3,46tm, 
M, =— 1,47tm. 


Das negative Vorzeichen zeigt, daB der Richtungssinn der Momentvektoren entgegengesetzt der 
in Abb. 3 bis 7 dargestellten Richtungssinn ist. | 


Es ist zu beachten, daB sich folgende Werte fiir die Biegemomente ergeben wiirden: 
ME, =— 0,61 tm, M, =— 2,01 tm, M,z = — 2,69 tm, M, = — 3,34tm, 
M,=—1,44tm, 
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wenn man bei der Aufstellung der Gleichungen simtliche Beiwerte Se (a). fy («) usw., die von der 
Torsionsenergie abhangig sind, vernachlassigen wiirde. 

Man ersieht aus dem Beispiel, da® der Unterschied fiir den ersten und zweiten Fall nirgends 
gréBer als 4% ist. Eine Vernachlassigung der geraden Beiwerte kann jedoch im allgemeinen nicht 
empfohlen werden, weil deren Einflu8 auf die Torsionsmomente viel gréBer als auf die Biege- 
momente ist. 

Nachdem die Biegemomente an den Stiitzen ermittelt sind, 
ist es leicht, mit Hilfe der fiir jeden Belastungsfall entwickelten 
Formeln die Stiitzentorsionsmomente sowie die Biege- und Tor- 
sionsmomente in einem beliebigen Querschnitt zu bestimmen, 
weil jedes Feld als ein selbstandiger Trager auf zwei Stiitzen 
(Abb. 1) betrachtet werden darf, der durch die angegebene duBere 
Belastung und die gefundenen Biegemomente an den Stiitzen be- 
lastet wird. 

So ergibt sich z. B. das Torsionsmoment, das die Stiitze 2 auf- 
zunehmen hat, auf Grund von (14), (7), (3) und (4) zu 


rf 1 & 1 
9 =e 1) Har? (tg a — a) — M, (— 2etg 2a) 


os & 
1 1 1 1 
as fe 20 za) M; e 2a z3) Sage aa 
Wenn man die geraden Funktionen vernachlassigt, dann erhalt 


man — 0,09tm. Das Biegemoment in einem beliebigen Quer- 
schnitt des Feldes 2—3 wird 

= : M,sin p 
M, = qr? (tg asin g — 1 + cos ¢) ate on 

+ M, (cos y — ctg 2 «sin gy) =2,0-3,0? (0,577 sin gp — 1 + cos—) 

— 3,180 sin p — 2,06 (cos p — 0,577 sin ~) = 8,40 sin p — 18,00 

+ 15,94 cosy, usw. 

4. Der festeingespannte Kreistriger als Sonderfall. Die er- 
haltenen Formeln fiir die Drehwinkel und die Zahlenwerte der 
Tabelle kénnen auch zur Berechnung von Kreistragern mit fest- 
eingespannten Enden verwendet werden. Es geniigt in diesem 
Falle, den Ausdruck fiir den Drehwinkel an der Einspannungs- 
stelle anzuschreiben und ihn gleich Null zu setzen. 

Als Beispiel betrachten wir einen festeingespannten Kreistrager 
mit 2 « —90°, der durch die Einzellast Pin der Mitte beansprucht AN, tO, 
wird (Abb. 9). 

Man schreibt den Ausdruck fiir den Drehwinkel in der Tangentialebene an der Einspannungs- 
stelle und setzt ihn gleich Null: 


Pr* Pr? Mr Mr 
ORF Ey -- 0,0706 age 0,7854 55 — 0,1488 te 0,5000 


P=5t g=2t/m 


Mr Mr 
wobei M das Biegemoment an der Stiitze ist. Daraus folgt 
0,2777 + 0,0706 t 
Wa ET aged 0, 285Ad 


Forint — ict. Beinn iat Af 0.2217+P rs 


Das Torsionsmoment an der Stiitze ergibt sich auf Grund von (14) (3) und (4) zu 


Pr/ 1 se 0,0115 + 0,0115 t 
PDE) ex i Mtga =— Pryyeeq0,20540" 
Bei t =1 ist m —— 0,0146 Pr. Der wirkliche Richtungssinn der Momentvektoren ist in Abb. 9 


dargestellt. 
Es ist zu bemerken, daB das hier entwickelte Verfahren auch zur Berechnung von Durchlauf- 


tragern verwendet werden kann, welche im Grundrif die in Abb. 10 gezeigte Form haben. 


(Eingegangen am 24, Oktober 1949.) 
Anschrift des Verfassers: Dozent Dr. Ing. W. Swida, Karlsruhe, Stuttgarter Str. 21. 


250 Grammel: Tafeln der verallgemeinerten Kreisfunktionen. Ingenieur-Archi} 


Tafeln der verallgemeinerten Kreisfunktionen 
Sin (4) v, Cos (4) v, Sin (6) v, Cos (6) v, Sin (8) v, Cos (8) v. 


Von R. Grammel. 


Ich bin mehrfach brieflich nach Tafeln der von mir vor kurzem ausfihrlich entwickelten very 
allgemeinerten Kreisfunktionen! Sin (n)v und Cos (n)v gefragt worden, welche da und dort fit 
die Lisung technischer Aufgaben verwendet werden, und lege daher hier solche Tafeln fitr dii 
wohl zunichst wichtigsten Falle n =4, 6 und 8 vor. Die ziemlich miihsame Berechnung diese} 
Tafeln ist nach einem von uns gemeinsam entwickelten Plan von meinem Mitarbeiter F. Jindr« 
sorgfaltig durchgefiihrt worden, und zwar mit sieben bis acht Stellen hinter dem Komma, wo) 
von nur die letzte und in wenigen Fallen die vorletzte unsicher war, so daB wir die Richtigkeit de: 
hier mitgetcilten sechs Stellen gewadhrleisten kénnen. 


0 OZ OF 36 48 10 
Abb. 1. 


Die Werte der Tafeln 1 bis 3 fiir ‘Sin (4)v, Cos(4)v, Sin(6)v, Cos (6)v, Sin(8)v und Cos (8) ¥ 
sind fast alle mit Hilfe der von mir seinerzeit angegebenen Reihen? dieser Funktionen fiir da: 
Argument v — 0,00 bis v= 1,00 gerechnet; nur fiir die Bereiche schlechterer Konvergen 
[v = 0,95 bis v = 1,00 bei Sin (6)v und Cos (6)v; sowie v = 0,90 bis v = 1,00 bei Sin(8)v une 
Cos (8) v] sind die Werte durch numerische Berechnung der Integrale fiir ArSin x und ArCos +. 
und Funktionsumkehr ermittelt worden. An sich hatte es geniigt, die Tafeln 1, 2 und 3 je bis 


zum Argument v =—-74(= 0,927 ...), v= 4 15 (= 0,964...) und v= 4-1 (= 0,978 ...) 


fortzusetzen, da mit den Tafelwerten von Cos (n)v und Sin (n)v fiir den Bereich v = 0 bis v= +m 


die beiden Funktionen fiir gerades n vollstandig bestimmt sind. Wir haben die Tafeln aber darithex 
hinaus (genauer gesagt iiber vy = 1 hinaus) in allerdings gréferen Schritten bis zum Argument 
v = 4,0 fortgesetzt, so daB sie mindestens eine halbe Periode iiberdecken. Selbstverstandlic 


sind diese weiteren Tafelwerte vermittels Spiegelung des Arguments an den Werten v = = 


1 
und > % und durch Interpolation gefunden worden. 


Weil es fiir manche Zwecke bequemer sein mag, das Argument in der Form 
= 
1 R. Grammel, Ing.-Arch. 16 (1948), S. 188. 
> Der Reihenkoeffizient 5”) (27) auf S. 193 ist zu berichtigen in: 


(n) _ eat 1 n 1 n)z n n 
oS ea 2) aq m — 3) (mn — 4) af +S (n— 3) af? a") +. al”) ah) + sigs tail, 
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Tafel 1. 
v | Sin (4) v Cos (4) v v | Sin(4) v Cos (4) v v Sin (4) v Cos (4) v 
0,00] 0.000000 | 1,000000 0,50 | 0.495388 | 0,984591 1,0 | 0,880880 | 0,794227 
01 | 010000 | 1,000000 51 504915 983 340 1 92 a0 
; 4.567 72 
02 020000 | 1,000000 52 514 405 982. 016 2 956 038 636 45 
03 030000 | 1,000000 53 523 855 980617 3 976 925 546 431 
04. 040000 | 0,999999 54 533 263 979 140 4 989 473 451211 
0,05 | 0.050000 | 0.999998 0.55 | 0,542628 | 0,977583 1,5 | 0,996085 | 0,353 243 
06 060.000 999 997 56 551 948 975 943 6 998 959 253 916 
07 070 000 999 994 57 561219 974.219 7 999 859 154.062 
08 080000 999 990 58 570 439 972 407 8 999998 | 0,054075 
09 089 999 999 984 59 579 608 970 506 9 999999 |—0,045 925 
0,10 | 0,099999 | 0,999975 0,60 | 0.588721 | 0,968512 2,0 | 0,999887 |—0,145 915 
11 109 998 999 963 61 597176 966 424. 1 999086 |—0,245 791 
12 119 996 999 948 2 606 772 964 239 2 996 432 |—0,345 185 
13 129 994 999 929 63 615 106 961 955 3 990202 |—0,443 308 
“4 139 992 999 904 64 624574 959570 4 978227 |—0,538 818 
0,15 | 0.149989 | 0,999873 0,65 | 0.633376 | 0,957082 2,5 | 0,958109 |—0,629 798 
16 159 984 999 836 66 642 107 954. 488 6 927572 |—0,713 887 
17 169 979 999 791 67 650 767 951 786 7 884919 |—0,788 618 
18 179972 999 738 68] 659351 948 975 8 829442 |—0,851 901 
19 189 963 999 674 69 667858 946 053 9 161635 |—0,902527 
0,20 | 0.199952 } 0.999600 0,70 | 0.676285 | 0,943017 3,0 | 0,683091 |—0,940458 
21 209 939 999514 71 684 629 939 866 1 596105 |-—0,966 818 
22 219 923 999 415 71 692 889 936 598 2 503155 |—0,983577 
23 229 904 999 301 73 701 060 933 213 3 406 463 |—0,993 105 
24. 239 881 999 171 74, 7109 142 929 107 4 307734 |—0,997 750 
0.25 | 0,249854 | 0,999024 0,75 | 0,717132 | 0,926080 3,5 208091 |—0,999531 
26 259 822 998 859 16 725 026 922 330 6 108147 |—0,999 966 
27 269 785 998673 17 732 823 918457 7 | 0,008149 |—1,000000 
28 279 142 998 465 78 740 521 914.459 8 |—0,091850 |—0,999 982 
29 289 693 998 235 79 748 116 910 335 9 |—0,191812 |—0,999661 
0,30 | 0,299636 | 0.997979 0,80 | 0,755608 | 0,906085 4,0 |—0,291534 |—0,998 189 
31 309572 997696 81 162 993 901 707 
32 319 498 997385 82 710 270 897201 Tafel La. 
33 329 415 997043 83 111437 892 566 
34 339 321 996 669 84. 184.492 887 803 : Ta Ta 
0,35 | 0,349215 | 0,996 261 0,85 | 0,791432 | 0,882 910 9 |Sin (4)@5*\Cos(4) 05 
36 359097 995 817 86 7198 256 877 888 
Sui 368 965 995 334 87 804 963 872 idea 0,00 {0,000 000 |1,000 000 | 1,00 
38 378 818 994 811 88 811550 867457 05 | 092703 |0,999982 | 95 
39 388 655 994 246 89 818017 862 047 10 | 185375 | 999705 | 90 
15 | 277862 | 998506 | 85 
0,40 | 0,398474 | 0,993637 0,90 | 0.824362 | 0,856509 20 | 369768 | 995293 | 80 
42] 418056 992 275 92 | 836681 845 050 30) | 0948433 | 9761973) |) «10 
43} 427815 991518 93 | 842652 839 129 35 | 632434 | 957354 | 65 
44] 437550 990 708 94 | 848.498 833 083 40 | 710451 | 929124 |) 60 
0,45 | 0.447262 | 0,989842 0,95 | 0,854216 | 0,826 912 45,3) SB 00S pia) 18s 
47 466 604 987933 97 865 268 814 200 iz = 
48 476 231 986 885 98 870 601 807662 Sin (4)o 
4.9 485 826 985 772 99 875 805 801 003 
7 
0,50 | 0,495388 | 0,984591 1,00 | 0,880880 | 0,794227 = = 1854075 


| 
| 
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Tafel 2. 
) Sin (6) v Cos (6)u 0) Sin (6)v Cos (6) v v Sin (6) v Cos (6) v 
0,00] 0,00000 1,000 000 0,50 | 0.499075 0,997408 1,0 0,910 046 0,869 442 
01 010 000 1,000 000 51 508 938 997083 1 951004 7199 025 
02 020 000 1,000 000 52 518 784. 996 724 2 976 300 715 379 
03 030 000 1,000 000 53 528 612 996 330 3 989 994 623147 
04 040 000 1,000 000 54. 538 419 995 898 4 996 427 526 276 
0,05 | 0,050 000 1,000 000 0,55 | 0,548 204 0,995 424 1,5.] 0,998 983 0,427311 
06 060 000 1,000 000 56 557965 994.906 6 999 794. 327573 
07 070000 1,000 000 57 567700 994 342 7 999 977 227617 
08 080000 1,000 000 58 577405 993 726 8 0,999 999 127621 
09 090 000 1,000 000 59 587079 993 057 9 1,000 000 0.027621 
0,10 | 0,100000 1,000 000 0,60 | 0,596 719 0,992 330 2,0 1,000000 |—0,072 379 
11 110000 1,000 000 61 606 323 991542 1 0,999996 |—0,172378 
12 120 000 1,000 000 62 615 887 990 690 2 999 932 |—0,272 365 
13 130 000 0,999 999 63 625 408 989 769 3 999556 |—0.372 261 
14 140 000 0.999 999 64 634 883 988 775 4 998154 |—0,471 756 
0.15 | 0,150 000 0,999 998 0,65 | 0,644309 0,987 704 2,5 0.994200 |—0,570011 
16 160 000 999 997 66 653 682 986 553 6 985009 |—0,665 206 
17 170000 999 996 67 662 999 985 316 7 966682 |—0,754155 
18 179 999 999 994 68 672255 983 989 8 934792 |—0,832 440 
19 189 999 999 992 69 681447 982 568 9 886025 |—0,895 642 
0,20 | 0,199 998 0,999 989 0,70 | 0.690570 0,981 049 3,0 0,819987 |—0,941392 | 
21 209 998 999 986 71 699 621 979 4.26 1 | 739537 |—0,970666 | 
22 219997 999 981 72 108 594 977695 2 649230 |—0,987111 | 
23 229 996 999 975 13 117486 975 851 3 553325 |—0,995158 | 
24 239 995 999 968 14 7126 292 973 891 4 454762 |—0,998520 | 
0,25 | 0,249 993 0,999 959 0,75 | 0,735 006 0,971 808 3,5 0.355158 |—0,999665 | 
26 259 990 999 949 16 143 625 969599 6 255234 |—0,999954 | 
27 269 987 999 935 Va 152144 967 259 q 155242 |—0,999 998 
28 279 984. 999 920 18 160557 964 783 3 0,055 243 |—1,000 000 
29 289 980 999901 719 168 860 962 168 9 |—0,044 757 |—1,000 000 
0,30 | 0.299974 0,999 879 0801 0,777048 0,959 408 4,0 |—0,144757 |—0,999998 | 
31 309 967 999 852 81 185 116 956 500 
32 319 959 999 821 82 793 059 953 441 Tafel 2a. 
33 329 949 999 785 83 800 872 950 225 
34, 339 937 999 743 84. 308551 946 850 Ts! Cos (6) 97S 
0,35 | 0,349 923 0,999 694 0,85 | 0.816091 0,943 312 3 |498 (Joa 
36 359 907 999 637 86 823 488 939 609 
37 369 887 999 573 87 830 738 935 738 0,00 0,000 000 |1,000 000 | 1,00 
39 389 837 999 414 89 844.777 927480 ‘ 
10 | 192761 |0,999 991 90 
0,40 | 0,399805 | 0,999318 0,90 | 0, 15 | 289123 | 999903 | 85 
we ta, Oe eeey eens ane e 20 | 385374 | 999453 80 
Al 409 769 999 209 91 858181 918 524 0,25 |0,481190 |0,997920 | 0,75 
42 419 726 999 087 92 864 636 913 780 30°} 575 744 | 993 835 710 
43 429677 998 948 93 870 922 908 859 35 | 667327 | 984 708 65 
AA 439 621 998 793 94 8717037 903 759 40 | 753031 | 967005 60 
0.45 | 0,449 556 0,998 619 0,95 | 0,882 980 0,898 481 45 | 828888 | 936749] 55 
46 459 483 998 425 96 888 747 893 024 0,50 10,890 899 |0,890 899 | 0,50 
AT 469 399 998 209 97 894 338 887391 5 
48 A719 304 997969 98 899 752 881582 Cos(6) 68! Sin(6) 5S 
49 | 489196 997 703 99 | 904988 875 598 eg pe ote e 
0,50 | 0.499075 0,997408 1,00 | 0,910 046 0,869 442 Ze = 1,927621 
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Tafel 3. 
a a ec ae ener ee 
v Sin (8) v Cos (8) v v Sin (8) v | Cos (8) v 
0,00 | 0,000000 | 1,000000 0,50 | 0,499810 | 0,999512 1,0 | 0,928220 | 0,904 736 
ol 010000 | 1,000000 51 509774 999 429 1 966 843 835 053 
02 020000 | 1,000 000 52 519730 999 333 2 987003 749 301 
03 030000 | 1,000000 53 529680 999 223 3 995 715 654 735 
04 040000 | 1,000 000 54] 539621 999098 4 998 847 556 422 
0,05 | 0,050000 | 1,000 000 0,55 | 0,549554 | 0,998 956 1,5 | 0,999 763 | 0,456 837 
06 060000 | 1,000000 56 559 475 998 795 6 999 967 356 912 
07] 070000 | 1,000000 57] 569385 998 612 7} 0,999 998 256 921 
08 | 080000 | 1,000000 58] 579281 998 406 8 | 1,000000 156 922 
09 | 090000 | 1,000 000 59 589 162 998 174 9 | 1,000000 | 0,056 922 
0,10 | 0,100000 | 1,000000 0.60 | 0.599026 | 0,997912 2,0 | 1,000000 |—0,043 078 
11 110000 | 1,000 000 61 608 871 997619 1] 1,000000 |—0,143078 
12 120000 | 1,000 000 62 618 695 997291 2} 0,999998 |—0,243 078 
13 130000 | 1,000 000 63 628 494 996 924 3 999976 |—0,343 072 
~ 14 140000 | 1,000 000 64 | 638267 996 515 4} 999814 |—0,443014 
0,15 | 0,150000 | 1,000 000 0,65 | 0.648010 | 0996060 2,5 | 0,999057 |—0,542 680 
16 160000 | 1,000 000 66 657721 995 554 6 996 380 |—0,641 270 
17 170000 | 1,000 000 67] 667395 994.993 7] 988733 |—0,736 590 
18 180000 | 1,000000 68 677029 994 373 8] 970586 |—0,823 969 
19 190000 | 1,000 000 69 686 619 993 687 9 934871 |—0,896 306 
0,20 | 0.200000 | 1,000 000 0,70 | 0696160 | 0,992932 3,0 | 0,877151 |—0,947499 
21 210000 | 1,000 000 71 105 648 992 100 1 7199636 }—0,977395 
22 220000 | 0,999999 ie 715 078 991187 2 109 280 |—0,991 760 
23 230000 | 0,999999 73 724 445 990 186 3 612649 |—0,997497 
24 240000 | 0,999 999 74 733 142 989 091 4} 513601 |—0,999 394 
0,25 | 0,250000 | 0,999998 0,75 | 0,742965 | 0,987894 3,5 | 0,413809 |—0,999 892 
26 259 999 999 997 16 752 106 986 589 6 313 840 |—0,999 988 
27] 269999 999 996 71 761159 985 169 7] 213843 |—0,999 999 
28 279 999 999 995 78 170117 | 983 626 8 113 843 |—1,000 000 
29 289 999 999 994 79 118 974 981953 9 | 0,013843 |—1,000 000 
0,30 | 0.299998 | 0,999 992 0,80 | 0,787721 | 0,980 143 4,0 |—0,086 157 |—1,000 000 
a1 309 997 999 989 81 196 351 978 187 
32 319997 999 986 82 804 857 916 078 Tafel 3a. 
33 329 995 999 982 83 813 231 973 Be 
34 339 994 999 978 84 821 464 9713 78) 06 (8) (78 
0,35 | 0.349992 | 0,999972 0,85 | 0.829548 | 0,968 758 Cem Din (Berg COst SOs, 
36 359 990 999 965 86 837476 965 962 
37] 369987 999 956 87 845 240 962977 9,99 }0,000000 {1,000 000 | 1,00 
38 379 984 999 946 88 852 832 959 7196 05 | 097846 |1.000000 | 95 
39 389 980 999 933 89 860 244 956 412 10 195 692 |1,000 000 90 
Fl ae 15 | 293537 |0,999993 | 85 
0,40 | 0.399975 | 0.999918 0,90 | 0,867470 | 0,952 821 20 | 391363 |0,999931 | 80 
Al 409 968 999 900 91 874 502 949018 0,25 |0,489075 |0,999 590 | 0,75 
42 419 960 999 879 92 881 334 944 997 30 | 586275 | 998245 | 70 
43 429 951 999 854 93 887 960 940 756 85 | 681755 | 994044 | 65 
im 44 439 940 999 824. 94. 894 375 936 291 AOE Pied peo O00 
0,45 | 0.449926 | 0,999 790 0,95 | 0,900574 | 0,931 600 45 | 853293 | 959593 | 55 
46 459 910 999 750 96 906 553 926682 0,50 |0,917004 |0,917004 | 0.50 
47| 469891 999 703 97} * 912310 921536 ra a 
A8 479 869 999 648 98 917841 916 162 Gos(8)0-5|Sin(8)o5] 2 
49 489 842 999 585 99 923 144 910 562 
av 
0,50 | 0.499810 | 0,999512 1,00 | 0,928220 | 0,904 736 > = 1,956 922 
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zu beniitzen, so sind in den Tafeln la bis 3a die Funktionswerte auch noch fiir 9 = 0,00 Mi 


y 


o = 1,00 dargestellt. (Vgl. Abb. 1, wo die Funktionen im Vergleich mit sinv und cosv (n = 2)) 
dargestellt sind.) Wegen der groBen Argumentschritte dieser Tafeln wird man im allgemeinen| 
nicht mehr mit linearer Interpolation auskommen, wenn man fiir Zwischenwerte auch nur drei. 
bis vier giltige Stellen haben will, sondern quadratisch interpolieren miissen. 

Um die berechneten Tafelwerte zu sichern, haben wir jede mégliche Nachpriifung vorgenom-| 
men: Bildung der ersten und héheren Differenzen; Uberwachung der Formel 


Cos"(n) v + Sin” (n) v= 1 


fiir jeden fiinften Tafelwert; Berechnung der Tafelwerte fiir v = 7% gemaB der Formel 


To. r, 1 \l/n 


polationsformel durchgefiihrt ?. 


1 Vel. R. Grammel, a. a. O. Formel (23), worin der Druckfehler: k? sn/dn statt k sn/dn zu berichtigen ist.) 
2 Vel. C. Runge u. H. Kénig, Vorlesungen iiber numerisches Rechnen, S. 108ff., Berlin 1924; oder Fry 
A. Willers, Methoden der praktischen Analysis, 5. 78ff., Berlin 1928. 


(Eingegangen am 1. Dezember 1949.) 
Anschrift des Verfassers: Prof. Dr. Richard Grammel, Stuttgart-N, Robert-Bosch-Strafe 101. 
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Die Berechnung der Auftriebsverteilung elastisch verdrehbarer 
Tragfliigel. 


Von J. Weissinger. 


1, Einleitung. Die Aufgabe, mittels des Prandilschen Traglinienmodells die Auftriebsvertei- 
lung eines starren Tragfliigels zu bestimmen, lat sich mathematisch bekanntlich auf zwei ver- 
schiedene Arten ausdriicken: erstens nach Prandtl! in Form einer Integrodifferentialgleichung 
und zweitens nach Trefftz? als potentialtheoretisches Randwertproblem. Zu diesen beiden For- 
mulierungen ist in neuerer Zeit noch eine dritte getreten: Ziller® erkannte, da® sich das Problem 
auch als Variationsaufgabe aussprechen laft. Zur Lisung dieser Aufgabe verwendet Ziller das 
Ritzsche Verfahren unter Benutzung eines trigonometrischen Ansatzes und gelangt so zu einem 
linearen Gleichungssystem fiir die gesuchten Fourierkoeffizienten der Auftriebsverteilung. Da 
man sich aber leicht iiberzeugen kann, dafB dieses Gleichungssystem mit dem von I. Lotz* auf- 
gestellten iibereinstimmt®, so ist das Zillersche Verfahren im Hinblick auf die numerische Durch- 
fiihrung mit dem Lotzschen identisch. Die Hauptarbeit bei diesem Verfahren macht bekanntlich 
nicht die Auflésung, sondern die Aufstellung des Gleichungssystems, da hierzu eine har- 
monische Analyse von zwei vorgegebenen Funktionen (im wesentlichen Anstellwinkel- und rezi- 
proke Tiefenverteilung) vorangehen muB. Aus diesem Grunde hat sich heute das vor einigen 
Jahren von Multhopp® gefundene Verfahren wohl ziemlich allgemein eingebiirgert, bei dem das 
lineare Gleichungssystem unmittelbar gegeben ist, so daf die einzige Arbeit in der durch Iteration 
leicht zu bewerkstelligenden Auflésung dieses Systems liegt. AuBerdem erhalt man als Ergebnis 
nicht die Fourierkoeffizienten, sondern unmittelbar die Werte der Auftriebsverteilung an be- 
stimmten Stellen der Spannweite. 


Im zweiten Teil seiner Arbeit hat Ziller gezeigt, daB sein Verfahren auch auf den elastisch 
verdrehbaren Tragfliigel, bei dem also die geometrische Anstellwinkelverteilung durch die Auf- 
triebsverteilung beeinfluBt wird, angewandt werden kann. Die numerische Durchfihrung aller- 
dings ist auSerordentlich miihsam, da zu den oben erwahnten noch zwei weitere harmonische 
Analysen hinzukommen, darunter eine nach zwei Variablen. AufSerdem hat das Verfahren den 
Nachteil, da8B es nicht aligemein anwendbar ist, sondern nur dann, wenn der Abstand des Neu- 
tralpunktes von der elastischen Achse lings der Spannweite konstant ist. Da andrerseits keine 
bessere Methode bekannt zu sein scheint’, soll im folgenden ein Verfahren geschildert werden, 
das gegeniiber dem Zillerschen drei Vorteile hat: 

1. Es ist unbeschrankt, also fiir beliebige ungepfeilte Fliigelformen und beliebige Lage der 
elastischen Achse anwendbar. 

2. Es fubt unmittelbar auf dem zurzeit woh] meist gebrauchten Verfahren von Multhopp und 
kann daher miihelos von jedem, der dieses beherrscht, benutzt werden. 

3. Der Rechenaufwand ist verhaltnismaBig gering, namlich nicht wesentlich gréfBer als der 
des gewéhnlichen Multhoppschen Verfahrens. 

Bei der letzten Angabe ist natiirlich (ebenso wie bei Ziller) vorausgesetzt, das die von der 
Statik zu liefernden Unterlagen, insbesondere also die TorsionseinfluBlinien des Fliigels, bekannt 
sind, wie wir uns iiberhaupt im folgenden auf die aerodynamische Seite der Aufgabe beschranken 
wollen. In Ziff. 3 wird das Verfahren allgemein entwickelt, in Ziff.4 der Rechnungsgang 


an einem praktischen Beispiel erlautert. 


1 ZL. Prandtl, Nachr. Ges. d. Wissensch., S. 451. Géttingen 1918, §. 451. 

2 E. Trefftz, Z. angew. Math. Mech. 1 (1921), 5S. 206. 

3 F. Ziller, Ing.-Arch. 11 (1940), 5. re 

4 J. Lotz, Z. Flugtechn. Motorl. 22 (1931), 5. 189. 

5 Bis auf einen konstanten Faktor sind namlich die Zillerschen B,, gleich den bei Lotz [Gleichung (4)] 
auftretenden Differenzen zweier Fourierkoeffizienten y,, und die C,, gleich den a, ,,. 

6 H. Multhopp, Luftf.-Forschg. 15 (1938), S. 153. 

7 Mir ist hie noel das Verfahren von G. Datwyler, Flugwebr u. -Technik Bd. 3 (1941), 5. 65, bekannt, 
das aber nur mit zwei Punkten der Halbspannweite arbeitet und daher im allgemeinen zu ungenau sein 
diirfte. 

gs 
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2. Bezeichnungen. Wir fiihren folgende Bezeichnungen ein: 


y Koordinate in Richtung der Spannweite (y = 0 = Fliigelmitte ; 


yor 2 = Fliigelenden); 


2 ‘ ; : 
i= ae dimensionslose Koordinate; 


0 = are cos 7 Winkelkoordinate ; 
b Fligelspannweite ; 
F Fliigelflache ; 

b? 
of is ¥F 

l(y) Fligeltiefe ; 
V Fluggeschwindigkeit ; 

o Luftdichte; 


= + V? Staudruck; 


Fliigelstreckung ; 


Gol¥) == Aty) Beiwert des értlichen Auftriebs A(y) (pro Langeneinheit der Spannweite); 


q !(y) 
cm (y) = ae Beiwert des értlichen Langsmoments (je Langeneinheit der Spannweite) bezogen H 
auf die elastische Achse im Profilschnitt y, schwanzlastig positiv; 
cm, (y)= aC Beiwert des ortlichen Langsmoments bei ¢,(y) = 0; 


A . ; ae 
6, == ae Beiwert des Gesamtauftriebs A ; 


a(y) Ortlicher geometrischer Anstellwinkel (gegen die Nullauftriebsrichtung des Profil- 
schnitts) im Bogenmab; 


dc 
Ci Profilkonstante E a | 
dasAa—>o 


I’(y) Zirkulation um einen Profilschnitt; 

pee tee oa) 
Nitya wap 

E(y) Entfernung des Profilneutralpunktes von der Profilnase dividiert durch I(y); 


é(y) Entfernung der elastischen Achse von der Profilnase im Profilschnitt dividiert 
durch I (y); 


O Torsionssteifigkeit ; 


dimensionslose Zirkulation; 


x ein die elastischen Higenschaften einer Strebe charakterisierender Parameter ; 


E(y, 9’) Einflu8funktion (TorsionseinfluBlinien) ; 


2 
EEA (y oy) E(y, y’) dimensionslose EinfluBfunktion; 


Bb(y) Verdrehwinkel ; 
Bb (y) Von M, (y) herrithrender Bestandteil des Verdrehwinkels. 


3. Darstellung des Verfahrens. Wir betrachten einen Fliigel von der Spannweite b und der 
Tiefenverteilung /(y). Seine elastischen Eigenschaften seien in Form von TorsionseinfluBlinien 
E(y, y’) bekannt. Bezeichnet dann M(y) irgendeine Momentverteilung mit der elastischen | 
Achse als Bezugspunkt, deren Entfernung von der Vorderkante im Profilschnitt an der Stelle ¥ 
durch €(y) /(y) gegeben sei, so wird an einer beliebigen Stelle y der folgende Verdrehwinkel her- | 


vorgerufen: 


b/2 


BOS] Ey) My) dy’. (1) 
Wird der Flugstaudruck mit q, die Riicklage des Neutralpunktes mit & (y) 1 (y) bezeichnet, und 
driickt man M durch den im Neutralpunkt angreifenden Auftrieb A und das bei A = 0 herr- | 
schende ,,Sturzflugmoment‘‘ Mj, sowie weiter A durch die dimensionslose Zirkulation y = I7/bV | 


( 
| 
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aus 


M(y) = A(y) Uy) [e(y) — &(y)] + Moly) = 298 y(y) Uy) Le (x) — §(y)] +4 emo(y) 2(y), (2) 
so geht (1) iiber in 


b/2 
Bly) = 2 9b f (ely) — Fy Uy) EU, V7) dy’ +89) (3) 
-wobei der Winkel 
: b/2 
BO (y) = 4} emo(y’) Ply’) Ely, y’) dy’ (4) 


von der unbekannten Auftriebsverteilung nicht abhangt und daher zunichst als bekannt an- 
gesehen werden soll. Fiihrt man schliefSlich noch die Variablen 77 und # mittels 


1= y= c08 9 (5) 


ein, so erhalt man (unter Beibehaltung der alten Funktionszeichen auch fiir die neuen Variablen) 
1 


B(y) = qb J [é(7’) — €(7’) U7’) Ey, 17’) y (7') dy’ +B (mn), (6) 


al 
ay fox o(77’) P(y') E(y, 1) dy . | (7) 
1 


ick 


Das in (6) auftretende Integral kann nun sofort nach der auch von Multhopp benutzten Qua- 
draturformel 


S faa = 5 aS fl) sin by (Oo= SG) (8) 


ausgewertet werden. Diese Gleichung ist exakt richtig, wenn f (7) gleich dem Produkt von ] 1—7? 
mit einem Polynom (2m —1)-ten Grades ist; sie liefert daher eine gute Naherung des exakten In- 
tegralwertes, wenn der Integrand sich als Produkt von ) 1—7/? mit einer stetigen Funktion dar- 
stellen 1aBt. Diese Voraussetzung ist aber fiir den Integranden von (6) erfiillt, da sie fiir y(7’) 
zutrifft und die iibrigen Faktoren in allen praktisch vorkommenden Fallen stetig sein werden. 
Bezeichnet man also die Funktionswerte im Punkte 7,= cos yze|(m--1) abkiirzend durch den 
Zeiger y und entsprechend den Wert E (n,, 7,) mit E,,, so darf man naherungsweise fiir (6) 
schreiben 


B= 2 ton Vn gO) (9) 

mit 
ay, = —— qb’1,.{@, — €,) Ey, sin 0, (10) 
Damit ist bereits alles Notige fiir die Anwendung des Multhoppschen Verfahrens vorhanden. 
Die elastische Verdrehbarkeit des Fliigels hat zur Folge, daB® in der Prandtlschen Integrodiffe- 
rentialgleichung der Anstellwinkel « um den Verdrehwinkel f vermehrt werden muf; d. h. statt 


des im starren Fall mit der Prandtlschen Gleichung aquivalenten Multhoppschen Gleichungs- 
systems 


by Y»= ay 4d) bon Yn (y=, 2pm os m) (11) 


i 


bekommt man jetzt 
bpp Oy Py On Vn > (12) 
n=1 
und wenn man hier den Ausdruck (9) einfiihrt, so ergibt sich 


be yy Oy Hop” $2 bias (of S35 152i iz wiegs:m) (13) 


mit 
b= by— Gy», Bes == bys + ayn . (14) 


| 
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Man erhalt also das neue Gleichungssystem (13) zur Bestimmung der y,, einfach dadurch, daf 
man in (11) den Anstellwinkel « durch « +6 und die Multhoppschen Koeffizienten 6 durch die 
b* nach (14) ersetzt, wobei die a nach (10) zu berechnen sind. Das System (13) kann auch nach 
derselben Iterationsmethode wie (11) aufgelést werden, vorausgesetzt, da der Verdrehwinkel 
nicht allzu groB wird. Der einzige formale Unterschied zwischen (11) und (13) besteht darin, daB 
eine Anzahl der b,, stets Null ist, wahrend die entsprechenden bj, nicht verschwinden miissen. 

Durch diesen Umstand ist eine kleine Modifikation in der numerischen Auflésung von (13) 
bedingt. Wahrend man namlich bei Multhopp zu Beginn der Iteration eine nullte Naherung le- 
diglich fiir die Werte mit geradem Zeiger n = 2n’ zu kennen braucht, ist bei (13) eine nullte 
Naherung fiir alle y, erforderlich. Fir gerade Zeiger bekommt man diese wie bei Multhopp durch 


Auflésung des Gleichungssystems mit = (m—1) Punkten, und indem man den Funktions- 


verlauf itber 7 an Hand dieser Werte aufzeichnet, kann man sich auch die fehlenden Punkte durch 
Ablesen beschaffen. 

Im Falle, daB® die Verteilung « + 8) symmetrisch oder antisymmetrisch (und natiirlich der 
Fliigel symmetrisch) ist, 1aBt sich das Gleichungssystem (13) ebenso wie das Multhoppsche 


weiter vereinfachen. Zu dem Zweck setze man 


Gin fir — oe : 
A(t) = 4 (15) 
Ayn | Oy, mo Tn fiir no, 
bzw. 
NN fae ee (16) | 
und bilde mit den Multhoppschen B,,, C,, die Koeffizienten 
By, == By, Sia) > (17) | 
bzw. | 
Chet G tas A Gels (18) | 
Dann lautet das Gleichungssystem im symmetrischen Falle | 
m-+1 
2 
BF y= a + HO +S BY yn ( =1,2,..., a) (19) 
n=1 
und im antisymmetrischen Falle 
m—1 


2 at | 
= mt HOLS Cre (r= 12. BEI). (20) 


Eine beliebige Verteilung lat sich in einen symmetrischen und einen antisymmetrischen Bestand- | 
teil aufspalten. | 
Die Kinflu8funktion E (7, 7’) hat im allgemeinen die Eigenschaft, daB sie fiir Wertepaare — 

n, 1’ mit entgegengesetztem Vorzeichen verschwindet. Dann vereinfacht sich (15) zu 
Att) =a . 


n 


(15a) 


Der ,,Nullverdrehwinkel* 6) (7) spielt fiir die Berechnung der Auftriebsverteilung dieselbe 
Rolle wie der Anstellwinkel a, ist allerdings im Unterschied zu diesem nicht unmittelbar gegeben, 
sondern muf} nach (7) erst berechnet werden. In den meisten Fallen wird sich das Integral (7) 
explizit ausrechnen lassen, sonst muB man sich die Werte 6) durch graphische oder numerische 
Integration beschaffen. 

SchlieBlich noch eine kurze Bemerkung iiber die Konvergenz des Iterationsverfahrens zur Lé- 
sung von (13). Wahrend sich fiir (11) die Konvergenz allgemein beweisen lat}, ist dies fiir (13) 
nicht méglich und auch nicht richtig ; denn wenn die a,, sehr groB sind, der Fliigel also sehr weich 
ist, so wird der Diagonalkoeffizient b} keineswegs mehr groB gegeniiber den b*, sein. Die An- 
wendbarkeit des Iterationsverfahrens priift man am besten mit Hilfe des bekannten »,Zeilen- 


1 Vgl. H. Multhopp, a. a. O. Abschn. VI. 
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kriteriums‘‘ 
Dr bs | < |b | Gis bh 25-2 4,2,.711) (21) 
p= 

nach. Entsprechendes gilt fiir die Gleichungssysteme (19) und (20). 


4. Ein Beispiel. Im folgenden soll das in Ziff. 3 geschilderte Verfahren an einem praktischen 
Beispiel erlautert werden. Der betrachtete Fliigel hat rechteckigen GrundriB und 14Bt sich als 
Torsionsréhre auffassen. Die Torsionsmomente in Fliigelmitte werden vom Rumpf aufgenom- 
men, dessen Breite vernachlissigt wird, so daB also die EinfluBfunktion E (y, y’) fiir Wertepaare 
y,y’ von entgegengesetztem Vorzeichen sowie fiir y— 0 oder y’=0 verschwindet. An der Stelle 


Ye 7) a (bzw. y = — yp) ist der Fliigel abgestrebt. 


Wenn man nur symmetrische oder antisymmetrische Verteilungen ins Auge faBt, aus denen 
sich ja eine beliebige Verteilung durch Uberlagerung gewinnen lat, so kann man sich von vorn- 
herein auf die Betrachtung einer Fliigelhalfte, etwa y 20, beschranken; und da E (y, y’) fiir 
y’ < O gleich Nullist, so darf man auch y’ >0 annehmen. Dann lebrt die Statik, daB E (y, y’) 
als Funktion einer der beiden Variablen (etwa y’) betrachtet, die Form eines doppelt geknickten 
Streckenzuges hat. Die Knickstellen liegen bei y’= y und y’= y,. Fiir die Rechnung empfiehlt 
sich die Einfiihrung der dimensionslosen Einflu8funktion 


E* (m9!) = 55 B (929) (22) 


wobei @ die Torsionssteifigkeit des Fliigels bedeutet. Indem wir schlieBlich noch eine von den 
geometrischen und elastischen Eigenschaften der Strebe abhingige Zah]l x einfiihren, erhalten 
wir auf Grund einer einfachen statischen Betrachtung die folgenden Definitionsgleichungen: 


a) fir 0 <9 <M: 


, i % , 
SU DEEE SY ESET B* (n,1!)=(1—“n) a, 
Br ea] a Yet =(l—f')a 
3.0 Sy SH S71’: = (1—~x)n, (23) 
b) fir 0 <m <7 
1.0 <7 <mx7: E*(n, 4’) = (1—x)y’, 
Ae Wh fA (OE = 7'— #N: 
35,0. on =n’ 3 ==. 1] 7 -# Noir 


Ist keine Strebe vorhanden, so ist x = 0 zu setzen. 


Da I(7) konstant ist und wir im folgenden auch nur konstantes ¢y(7)) betrachten wollen, so 
lauft die Berechnung von B® (7) nach (7) auf die Integration von E* (7, 7’) hinaus, die sich ele- 
mentar ausfihren lat. Das Ergebnis lautet 


1 
a) fir 0 <n <m: f B+ na) anf =n "$1 %} | (24) 


1 
‘ ’ ’ “x — 1° 
b) fir 0 Sm <q: J E*(y, 1!) dy! = hn — # M0 
Damit sind alle fiir die Berechnung der Auftriebsverteilung notwendigen statischen Unterlagen 


bereitgestellt. 
Der der Rechnung zugrundegelegte Fliigel ist durch folgende Zahlenangaben gekennzeichnet : 


b = 14,50 m, 1=1,65m (also F= 23,42 m?; A= 8,98), 
Coa .6,0, é= 0,334, & = 0,235, Ng = 0,608 , 
O = 76300 kg m?, % = 0,1023. 


Verwindung ist nicht vorhanden. Der Flugstaudruck sei g= 565 kg m*%. Das Multhoppsche 
Verfahren wird mit m= 15 angesetzt. Seat 

Zunichst hat man nach (23) die Werte E¥, zu berechnen. Das Ergebnis ist in Tabelle 1 ent- 
halten. Man erkennt, daB es sich nicht um 64, sondern nur um 13 auszurechnende Werte handelt. 


262 Weissinger: Die Berechnung der Auftriebsverteilung elastisch verdrehbarer Tragfliigel. — Ingenieur-Archiv 
pid ita teicher: tactic tb ied tats x sles BLT SPR BR 


und der Auftriebsgradient d¢,/da des elastischen Fliigels ist durch 
dca 


= (h) 
Alige a (30) 
gegeben. 
Im vorliegenden Beispiel war fiir cm der Wert —0,03 vorgeschrieben und die Auftriebsvertei- | 
lung sollte (bei g = 565 kg/m?) zu den Beiwerten ¢, = 0 und ¢, = —0,197 berechnet werden. Das 
Ergebnis, das man auf Grund von (28) mit (0) = 1,7l und e) = 5,51 erhialt, istjin Abb. 4 und 5 


zusehen. Zum Vergleich sind auch die mit m= 7 errechneten Werte eingezeichnet, die sehr gut | 


mit der genaueren Rechnung iibereinstimmen, so daf man sich bereits mit ihnen hatte begniigen | 
kénnen. SchlieBlich ist das nach dem Zillerschen Verfahren gewonnene Ergebnis eingetragen?. | 


Die Reihen wurden nach dem Koeffizienten mit dem Zeiger 11 abgebrochen, also ein System von 


6 Gleichungen aufgelist. Auch hier ist die Ubereinstimmung sehr gut. Kleine Unterschiede | 


weisen beide Verfahren hinsichtlich der elastischen Verwindung auf (Abb. 6), die bei Ziller nicht, 


wie es sein muB, fiir 7 — 0 verschwindet und im iibrigen eine leichte Welligkeit zeigt. Der Grund | 


dafiir ist wohl in der Approximation von f (7) durch eine abgebrochene Fourierreihe zu suchen. 


5. Zusammenfassung. Die Arbeit untersucht die Auftriebsverteilung des elastisch verdreh- | 
baren Tragfliigels, und zwar den aerodynamischen Teil der Aufgabe, setzt also die elastischen | 


Eigenschaften des Fliigels als bekannt voraus. Zur Lésung der Aufgabe hatte Ziller ein Verfahren 


angegeben, das einen betrachtlichen Zeitaufwand erfordert und nur auf Fliigel mit konstantem | 
Abstand zwischen elastischer Achse und Neutralpunktslinie anwendbar ist. Das vorstehend be- | 


schriebene Verfahren vermeidet diese Nachteile und schlieSt sich auBerdem unmittelbar an das 
beim starren Fliigel allgemein gebrauchliche Multhoppsche Verfahren an. Der Rechnungsgang 
wird an einem Beispiel erlautert, die Ergebnisse stimmen mit den nach Ziller gewonnenen 
iiberein. 

Bei unstetigem Anstellwinkelverlauf wird das Verfahren in der beschriebenen Form (ebenso 
wie das Zillersche) zu langsam konvergieren, man mu dann vielmehr mit der von Multhopp an- 
gegebenen Abspaltung arbeiten. Auf die Wiedergabe des zugehérigen Formelapparates wurde 


verzichtet. Desgleichen wurden Fliigel mit nicht gerader 1/4-Linie (z. B. Pfeilfliigel) nicht be- 


trachtet. Sie lassen sich ganz analog behandeln, wenn man nicht wie hier die Prandtlsche Trag- 
liniengleichung, sondern eine vom Verfasser angegebene allgemeinere Gleichung® zugrundelegt. 


1 Die Rechnungen wurden von Herrn W. Richter ausgefiihrt. 
2 J. Weissinger, Math. Nachr. 2 (1949), S. 45. 


(Eingegangen am 9. Dezember 1949.) 
Anschrift des Verfassers: Dr. J. Weissinger, Hamburg, Harvestehuderweg 10, 
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Uber eine Anwendung der Fouriertransformation 
in der Elastizitatstheorie*. 


Von H. Jung. 


1. Einleitung. Die Berechnung der Spannungen und Formanderungen kann beim ebenen 
Spannungszustand sowie beim ebenen Verzerrungszustand auf die Bipotentialgleichung 


oF 4 oF 
AAF (x, y) =>] 2 52098 Gree (1) 


zuriickgefiihrt Werden, wobei F(x, y) die Airysche Spannungsfunktion ist. 


Im folgenden wird in Ziff. 2 gezeigt, daB sich durch die Anwendung der Fouriertransformation 
die Airysche Spannungsfunktion direkt aus den Randbedingungen berechnen 1a$t. Damit wer- 
den in Ziff. 3 die bekannten Lésungen fiir die Halbscheibe und den unendlich langen Streifen 
mit Normalbelastung auf neuem und einfacherem Weg hergeleitet. AuSerdem laBt sich so, wohl 
erstmalig, die Lésung fiir den unendlich langen Streifen mit Tangentialbelastung angeben (Ziff. 4). 
Wie daraus auf einfache Weise die Lésungen weiterer wichtiger Probleme gefunden werden kén- 
nen, Wird in Ziff. 5 am Beispiel des durch eine Einzelkraft im Innern belasteten Streifens gezeigt. 


2. Transformation der Differentialgleichung und der Randbedingungen. Liegt ein ebener Span- 
nungszustand vor, so erhalt man in bekannter Weise! die Spannungen aus der Airyschen Span- 
nungsfunktion F(x, y) zu 

_ er er _ eF ; 

On ayF? Cys 5 ae deme Y's >The (2) 

Die Spannungsfunktion ist ein partikulares Integral der Differentialgleichung (1), das die Rand- 
bedingungen befriedigt. 

Multipliziert man (1) und (2) mit e—‘4* und setzt nun voraus, dafs die Integrale iiber x von 
— © bis + © existieren, so ergibt sich 


400 + 00 gay 
: oF : oF : 
CE emttede 42 f ee, e—i4xdx + | He sde—o, (la) 
+ co ={00 
2 
[osettax = [PF eed 
+ co 4 09 
2 
[eset dx = ota da, (2a) 


Unter geeigneten Voraussetzungen iiber F(x, y) gilt: 


1. die Differentiation nach y ist mit der Integration nach x vertauschbar, 
2. bei den partiellen Integrationen nach x verschwinden die integralfreien Teile. 


Dann erhalt man aus (la) und (2a) 


++ 00 ++ 00 + 00 
y t oe i ot —ihx aoe 
7h | Fee nertae— 2k [ Fe nerves +35 [ Fen Ae x= 0, (3) 


* Meinem Lehrer, Herrn Prof. Dr. R. Grammel zum 60. Geburtstag gewidmet. 
1 C. B. Biezeno u. R. Grammel, Technische Dynamik, S. 140. Berlin 1939. 
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und der Auftriebsgradient d¢,/da des elastischen Fliigels ist durch 
dca 
da 


gegeben. 
Im vorliegenden Beispiel war fiir cy der Wert —0,03 vorgeschrieben und die Auftriebsvertei- 


a4 Gd (30) 


lung sollte (bei g = 565 kg/m?) zu den Beiwerten ¢, = 0 und ¢, = —0,197 berechnet werden. Das | 


Ergebnis, das man auf Grund von (28) mit ¢\°) = 1,71 und é) — 5,51 erhalt, istiin Abb. 4 und 5 


zu sehen. Zum Vergleich sind auch die mit m= 7 errechneten Werte eingezeichnet, die sehr gut | 
mit der genaueren Rechnung iibereinstimmen, so da8 man sich bereits mit ihnen hatte begniigen _ 


kénnen. SchlieBlich ist das nach dem Zillerschen Verfahren gewonnene Ergebnis eingetragen?. 


Die Reihen wurden nach dem Koeffizienten mit dem Zeiger 11 abgebrochen, also ein System von | 
6 Gleichungen aufgelést. Auch hier ist die Ubereinstimmung sehr gut. Kleine Unterschiede 


weisen beide Verfahren hinsichtlich der elastischen Verwindung auf (Abb. 6), die bei Ziller nicht, 


wie es sein muB, fiir 7 — 0 verschwindet und im iibrigen eine leichte Welligkeit zeigt. Der Grund | 


dafiir ist woh] in der Approximation von f (7) durch eine abgebrochene Fourierreihe zu suchen. 
5. Zusammenfassung. Die Arbeit untersucht die Auftriebsverteilung des elastisch verdreh- 


baren Tragfliigels, und zwar den aerodynamischen Teil der Aufgabe, setzt also die elastischen | 


Eigenschaften des Fliigels als bekannt voraus. Zur Lisung der Aufgabe hatte Ziller ein Verfahren | 


angegeben, das einen betrachtlichen Zeitaufwand erfordert und nur auf Fliigel mit konstantem | 
Abstand zwischen elastischer Achse und Neutralpunktslinie anwendbar ist. Das vorstehend be- | 
schriebene Verfahren vermeidet diese Nachteile und schlieBt sich auSerdem unmittelbar an das | 
beim starren Fliigel allgemein gebrauchliche Multhoppsche Verfahren an. Der Rechnungsgang | 
wird an einem Beispiel erlautert, die Ergebnisse stimmen mit den nach Ziller gewonnenen | 


iiberein. 


Bei unstetigem Anstellwinkelverlauf wird das Verfahren in der beschriebenen Form (ebenso | 
wie das Zillersche) zu langsam konvergieren, man mu dann vielmehr mit der von Multhopp an- 


gegebenen Abspaltung arbeiten. Auf die Wiedergabe des zugehérigen Formelapparates wurde 
verzichtet. Desgleichen wurden Fliigel mit nicht gerader 1/4-Linie (z. B. Pfeilfliigel) nicht be- 
trachtet. Sie lassen sich ganz analog behandeln, wenn man nicht wie hier die Prandilsche Trag- 
liniengleichung, sondern eine vom Verfasser angegebene allgemeinere Gleichung*? zugrundelegt. 


1 Die Rechnungen wurden von Herrn W. Richter ausgefihrt. 
2 J. Weissinger, Math. Nachr. 2 (1949), S. 45. 


(Eingegangen am 9. Dezember 1949.) 
Anschrift des Verfassers: Dr. J. Weissinger, Hamburg, Harvestehuderweg 10. 
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Uber eine Anwendung der Fouriertransformation 
in der Elastizitaitstheorie *. 


Von H. Jung. 


1, Einleitung. Die Berechnung der Spannungen und Formanderungen kann beim ebenen 
Spannungszustand sowie beim ebenen Verzerrungszustand auf die Bipotentialgleichung 
or oF or 
AA F (x, y) =u | 2 aay ay 
zuriickgefiihrt Werden, wWobei F(x, y) die Airysche Spannungsfunktion ist. 

Im folgenden wird in Ziff. 2 gezeigt, daB sich durch die Anwendung der Fouriertransformation 
die Airysche Spannungsfunktion direkt aus den Randbedingungen berechnen la8t. Damit wer- 
den in Ziff. 3 die bekannten Lésungen fiir die Halbscheibe und den unendlich langen Streifen 
mit Normalbelastung auf neuem und einfacherem Weg hergeleitet. AuBerdem laft sich so, wohl 
erstmalig, die Lésung fiir den unendlich langen Streifen mit Tangentialbelastung angeben (Ziff. 4). 
Wie daraus auf einfache Weise die Lésungen weiterer wichtiger Probleme gefunden werden kén- 
nen, Wird in Ziff. 5 am Beispiel des durch eine Einzelkraft im Innern belasteten Streifens gezeigt. 


= 0 (1) 


2. Transformation der Differentialgleichung und der Randbedingungen. Liegt ein ebener Span- 
- nungszustand vor, so erhalt man in bekannter Weise! die Spannungen aus der Airyschen Span- 
nungsfunktion F(x, y) zu 

oF Ege Pee 9 
Oy?’ Sia eS AES oy” (2) 
Die Spannungsfunktion ist ein partikuladres Integral der Differentialgleichung (1), das die Rand- 
bedingungen befriedigt. 


Multipliziert man (1) und (2) mit e—‘4* und setzt nun voraus, da die Integrale iiber x von 
— © bis + © existieren, so ergibt sich 


40 7 Cae 
oF —tsay 42 f OF etinde + | OR tiedn=0, (la) 
es +00 
[osetia [FP e—eae, 
i ae 
[oyetednm [PF e—teae, (2a) 
ee <i, 
[xe-edx=— ape tae. 


Unter geeigneten Voraussetzungen iiber F(x, y) gilt: 
1. die Differentiation nach y ist mit der Integration nach x vertauschbar, 
2. bei den partiellen Integrationen nach x verschwinden die integralfreien Teile. 


Dann erhalt man aus (la) und (2a) 


+ oc + co + © 
} t 0 i Ax ot —ihx == 
A! [ Fl nena 28 [ Fe ne# dx T3y F(x,y)e—*44dx== 0,.. (3) 


* Meinem Lehrer, Herrn Prof. Dr. R. Grammel zum 60. Geburtstag gewidmet. 
1 C. B. Biezeno u. R. Grammel, Technische Dynamik, 5. 140. Berlin 1939. 
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oo + co +00 sta 02, 


+ 
[oe-teae=% [ela eds, [ ope Orde=— 7 | rey) e—i4x dx, 
——) Fo) —oo — co —— oo (4) if 
--+ co -++ co 


[ roctnda = we [ Planer ax ; 
y 
Fiihrt man nun in (3) ind (4) die Fouriertransformation? 


+ 0 +00 


F(x, yye-!#*dx=f(a,y), = | f(ayett* d= F(my), 
J 4 , al (5) 


F(x, y) > f(y)» F(A, y) > F(x, y) 


in Zeichen 


ein, so erhalt man 


{ ‘4 
AA F(x,y) = 0-—48f(d.y) — 28 FEM 4 SLED — 0, (3a) 
o, > SO, Ge eA oi 7a ; (4a) | 


Durch die Transformation (5) wird die Funktion F(x, y) vom Oberreich (x, y) auf den Unter- | 
bereich (A, y) abgebildet. 

Sind auf Randern y-konst. die Spannungen vorgegeben und lassen sich diese Spannungen in | 
den Unterbereich transformieren, so ist durch (3a) und (4a) ein einfaches Randwertproblem ge- 
geben. 

Die allgemeine Lésung von (3a) ist bekannt und lautet: 


f(A, ¥) = [40(4) +y Ar A)le— *¥ +-[ Bo (A) + y By(A)] e’¥. (6) 
Die Funktionen A; (A) und B;(A) lassen sich aus den in den Unterbereich transformierten Rand- 
bedingungen ermitteln. 
Durch Riicktransformation (5) erhalt man aus (6) die Spannungsfunktion fiir die auf Randern 
y-konst. belastete Scheibe (— © < x < + 0) 


+ co 
if . 
NES aS ee fae + y 4,(A)] e~*¥ + [By (A) + y By (A)je*7} ef" d2 . (7) 
Aus (7) lassen sich rein formal die Spannungen angeben. Es ist 
+ 00 
on= 5x | [2fU, x) —24(Aye—47 — ByeX)jet**aa, 
2 ; a 
v1=— za | 22/0. 9)e*aa, (8) 
ae 
t= [PLA + yA Ale“ +2 (BA) Ly B Wey 


A A,e—*7 1-1 Bi ett} ete dA. 


Da einfach zusammenhangende Bereiche vorliegen (— 0c <x <-++c), kénnen die Verschie- 
bungen wu und vin der x- und y-Richtung aus den Beziehungen? 


Owinrd 1 
Oa (>. m 05] i | 
Co 2a 1 Or 
7 Be” | 
1 W. Magnus u. F. Oberhettinger, Formeln und Satze fiir die speziellen Funktionen der Mathematik und 


Physik, S. 160, 2. Aufl. Berlin 1948. 
2 Siehe FuBnote 2 auf S. 263. 


bestimmt werden. 
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Geht man mit (8) in (9) ein und integriert, so erhalt man die Verschiebungen 


+ co 


Wa sez | oS STG. ye [4,(@)e-*y —B,()e4] ed). 


(10) 


+ co 
1 
Sear fi (m— 1) [(4,(a) + y 4, (A))] Ae-42 — (_By(A) + y B,(A)) Ae4y 
— (m +1) [A,(A)e—*y + B, (A) e*7]} ef ** dd. 

Sind auf Randern y-konst. die Verschiebungen vorgeschrieben, so kénnen die Funktionen A; (A) 
und B; (A) aus (10) durch Umkehrung der Fouriertransformation bestimmt werden. 

Die durch (7) gegebene Spannungsfunktion eignet sich besonders zur Behandlung von Be- 
lastungsfallen der Halbscheibe. Fiir den unendlich langen Streifen nimmt man die Lésung von 
(3a) in der Form an 


F(A, ¥) = a (A) Sin Ay + a,(A) Cof Ay + y [ag (A) Sindy + ay (A) of Ay] « (6a) 
Aus (6a) ergeben sich dann die Spannungsfunktion und die Spannungen 


+ co 


F(x y)= gh [ la(a) Sin ay+ay(@) Gof ay + as(d)y Sindy +ala)y Cojaylet*da, (Ta) 


—c 
+00 


== = [A2f(A, y) + 2A(a;(A) CojAy + a,(A) Sindy) e4*da 7 


— oo 


oy = — 2 [ara. y) end), 


+00 
ee +a, (A) A Sindy +a, (A)A of Ay] ef4=dd. 


Fiir die Verschiebungen erhalt man entsprechende Ausdriicke. Bei den einzelnen Belastungs- 
fallen ist jeweils zu zeigen, daB die in (7), (8) und (10) auftretenden Integrale gleichférmig in x 
und y konvergieren, um die angegebene Herleitung der Lésung zu sichern. Der Vorteil dieser Me- 
thode gegeniiber der einseitigen Laplacetransformation ist der, dai durch die Fouriertransfor- 
mation eine eindeutige Riicktransformation gegeben ist. 

3. Die Halbscheibe und der unendlich lange Streifen unter Normalbelastung. Die elastische, 
isotrope Halbscheibe werde an ihrem Oberrand y = 0 durch eine Normalbelastung beansprucht 
(Abb. 1). Die Randbedingungen fiir diesen Belastungsfall sind 


sr) Oy = yy (x) , (11) 6,-ys (x) 
i UO) 8 == 0, (12) 


Transformiert man die Belastung y,(x) in den Unter- 
bereich, so gilt 


8 


Abb. 1. Halbscheibe mit Normalbelastung. 


“feo 
J@yede=—Pe). 08) 


— oo 


Mit (13) werden die Randwerte im Unterbereich gema (4a) 
d 
VRS 0) ii Sir 8) — 0. (14) 


Um handliche Integrale fiir die Spannungen zu bekommen, werde die Unterfunktion fiir die 
Spannungsfunktion wie folgt angesetzt : 


SOY) =f Gay -hU ¥) 
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mit 
_fA(Aje*7 +A (A)ye—4+ firA>0, 
fils ¥) =) 9° ; fir A< 0, 


f(A; y) = | By (A) e4¥ + B,(A) yey fiir A =< () . 


Aus den Randbedingungen (14) erhalt man 


Av = 9A); Ay AGA) (15); 
By= 9,(4), By =—Aq,(A). 


Geht man mit (15) und (13) in (8) ein, so werden die Spannungen 


co B 
os ~ ff fox 1) i(6) e204 46-2) de dd 
0 «a 


—o ff 


+ | [ayy eereae—azar), 


co B 
y= sal] | Ay +I yQe-20 +16») de aa | 
7A ME (16) | 
0 a | 


| 
Ait | 


~| far- 1) p,() eA —#E—2) teal 
0 a 


~ B 
r=—34| Ay yy (é)e—4 +#E—) dédA 
0 4a 


5/7) | 
— | farmseyer seve azaal | 
0 «a 


| 


In (16) vertauscht man die Integrationsfolge und erhalt die schon von Féppl abgeleiteten In- | 
tegrale1 


B 
ae : [iy ee pilé)dé, 
= 
i mile: Saar las, (16a), 
2 | 
peer ae (E — x)y? 


my bet e— ap Mls) ae 


Wird die Halbebene durch eine Momentverteilung belastet, so ist bekanntlich 


d M(é) 
9 (€) = — d& s 


Damit wird, wie durch partielle Integration einfach gezeigt werden kann, 


+00 
pl) = f M(g)e-H8 de, 


woraus sich nach (15) und (8) die Spannungen ergeben. 


1 L. Féppl, Drang und Zwang, Bd. 3, S. 26. Miinchen 1948. 
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Ein unendlich langer Streifen von der Héhe 26 werde durch Normalkrifte belastet. Die 
Randbedingungen (Abb. 2) hierfiir sind 


Vo=) Osea w(x), T= 0", (18) 
y= —0:6,—%,(%), .7=0. (19) 
Transformiert man (18) und (19) in den Unterbereich, so erhalt man mit 
+ 0c 
vl ; 
A ama ey | (i (E) — s(E)) e—*4F de, 
+2 
1 ; = = 
V1 +P3= — ra KG + y5(f)) ede, ee ENE, 
Abb. 2. Streifen mit Normalbelastung. 


— © 


die Funktionen 
Cof Ad +16 Sindd 
wi 6) Or Ps) ein ead — 04d... 


SinAd + A6 Cojrzo 


a,(A) = (P1 +9) Gin 2ad-p2a0. 7° 
A Sin Ad (21) 
a3(A) = — (~, + @s) Gin2A0+ 240” 
ACojAd 
a4(2) = — (Y1 — 9) Gino se 
Geht man mit (21) in (8a) ein, so sind die Verschiebungen und Spannungent 
+0 } 
— al (1 — s) Ks + (G1 +s) KylAet** da 
fe Sin Ad CofA 
i Coj Ad Sindy ut DJ) AY ihe 
~ ZE | laa) Sn 21s — 215 + Pt Ps) sie oe 
+c 
1 idx 
pa ahy | 9) Ke HOt 9s) Ke ae a2 
NF CofA Sin 2d Sind (22) 
il Coj Ad Cojay im uLAy idx 
—4/ a4 Gin 245— 218+ (Pr + %) ee lee ds 
+00 
1 tax 
el £ | ao) K+ +0) K,] Wei** di, 
) + 
! 1 idx 
| o=— 2, | l—oo K, + (y+ 93) Ki] Aei** dd, 
+ co 
r=— 35 [tma—o) Ks+(~it Gs) Ke] ei** dd, | 


Eco 


wobei zur Abkiirzung gesetzt ist 
_ (46 Gin 45 —GCof 16) Sindy —Ay Coj Ad Cof ry 


ue Gin 246 — 2/6 ; 
K, = Ud Goj45 = SinAd) Cofhy—Ay Gindd Sindy 
_ Gin 226+ 226 ; 


1 F. Seewald, Abhandlungen aus dem Aerodynmaischen Institut d. Techn. Hochschule Aachen 5. 11. 
Berlin 1927. 
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(26 Gin Ad + Cof Ad) Sry —Ay Coj Ad CojA 


Uae Gin 216 — 2/6 y 
Ke (26 oj Ad 4- Sin Ad) Coj Ay —/Ay Sin 16 ae 
te Gin 246+ 226 
K 16 Sin Ad CojfAy —Ay Coj Ard Sindy 
aye Gin 210—2106 : 
K, Ad EojA6 Sindy —Ay Sindd Co} 


Gin 226 + 2106 


Damit die Integranten in (22) fiir A > 0 endliche Grenzwerte besitzen, miissen bei der Potenz- 
reihen- ce ua nach 4 die Koeffizienten von A—” (n > 0) verschwinden. Dies ergibt 


[Wy ()eE=0, f(wl)—wOEde=0, (23) 


+ co 


+ co 
f (vil)— pl) )24dE=0,  f (yl) — ya (€)) dE = 0. (24) 
Die Integrale (23) sind die bekannten statischen Gleichgewichtsbedingungen. Werden die Glei- 
chungen (23) und (24) durch die Belastung erfiillt, so konvergieren die Integrale (22) gleichmaBig 
in x und y, so daB durch (22) die Lésung des Randwertproblems gegeben ist. 


4, Der unendlich lange Streifen mit Tangentialbelastung. Wirken auf den Randern y= +0 | 
eines unendlich langen Streifens tangential angreifende Krafte, so lauten die Randbedingungen | 


(Abb. 3) 


Y= ONO 5 ofa Ve ST =o, (2), (25) | 
VO. Oy = Oe t= y, (x) (26) 
Damit erhalt man in entsprechender Weise wie in Ziff. 3 die Spannungsfunktion. Setzt man 
Say 
L 
2+ = —;/ [p2(€) +yalEle— i dé , 
DL 56%) 
T= Yy (#) Oye u 
: i oy walle "8 ag 


Poel tea [we (€) — 


Abb. 3. Streifen mit Tangentialbelastung. 


so werden die Verschiebungen und Spannungen 
++ co 


i . 
Uae |: Pa) Ls + (P2— Yq) Lylaet** dd 
+ oo 
stg ke Gin2d Sindy CofASCofAy | sae 
ai G «) m2 — 240 r (Pe a) Gm2Ad 4 210 e ahe 
bes 
u . 
ae Sacer aG [Ps + Pa) Ls ar (@ 5 — Ya) Lelet** da 
nae CinAdGojry | Goj4d Sindy | sas 
mt | lr 4) ©in 2A6—216 (P2 Pa) Gin 226 +216 Ce dj, (27) 


8 


38 


\(P2 + Pq) Ly + (Po Sait Py) Le] Aet4* di, 


Q 

| 
4 oe 
——— 


+ co 
1 ; : 
Sige ee [2 + a) Ls + (P2— Pq) Ly A? et ** d2., 
-+ 00 


| =: 


[G2 + Ga) Ls + (P2— Gy) Le Aet** di. 


i) 
a 


8 
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Hierbei sind die Abkiirzungen verwendet 
feel Sin 1d Sindy +Ay Sin1d Gojdy —16Coj16 Sindy 
Gin 2160 —2106 : 
eee Seay onde Oe As Sindy —A6 Sin 1d Cofiy 
Gin 200 +246 
pA y¥ Sin Ad Coj Ay — 4 CofAd Sindy 
Sin 216 — 2/16 
1, = ySei2e Sindy — 6 SinsAd Cojay 


Gin 246+ 2/0 : 
Lee Sin Ad CofAy +Ay Sin Ad Sindy —16 Coj rd Cojay 
Gin 216 —2A6 
pt, — S46 Sindy +AyCoiAd Coj Ay —AS Sins Sindy 
¥ Gin 205+ 2106 


Endliche Grenzwerte der Integranden in (27) fiir A > 0 erhalt man (vgl. Ziff. 3), wenn die Be- 
lastung den Bedingungen geniigt 


boo 
Sta @+mOles=0, 


oh () +y(le ae —0, | (28) 
a 


Sy Ey w,(e)e7 de — 0, Abb, 40 Streiien mite Bivvollase san tamern: 


Die Integrale (27) sind damit gleichmaBig in x und y konvergent, so da® durch (27) die Lésung 
des Randwertproblems gegeben ist. 

Fiir die durch tangentiale am Rand y = 0 beanspruchte Halbscheibe wird die Lésung ganz 
entsprechend wie in Ziff 3 durchgefithrt. Daher unterdriicken wir die Rechnung im einzelnen. 

Sind an einem unendlich langen Streifen allgemeine Belastungen auf den Randern y = +-6 
und y — — 6 vorgegeben, so lassen sich mit der hier beschriebenen Methode die Spannungen 
und Formanderungen angeben. An einem speziellen Beispiel soll in Ziff. 5 der Rechnungsgang 
gezeigt werden. 


5. Der unendlich lange Streifen mit einer Einzellast im Innern. Ein unendlich langer Streifen 
werde durch eine Einzellast im Innern beansprucht und liege auf dem Rand y = 6 in zwei Punk- 
ten auf (Abb. 4). Mit den in Ziff. 3 und 4 angegebenen Gleichungen lassen sich fiir diesen Be- 
lastungsfall die Spannungen angeben. 

Fiir eine Einzelkraft zunachst, die im Koordinatenursprung der vollen Scheibe angreift, wird 
mit Abb. 5 die Spannungsfunktion? 


+ co 
see ola Ee. 1 Lisa \ Sy ee y 
r= = qi E [ Ala eb **sinayda. (29) p 
Aus (29) erhalt man die zugehérigen Spannungen J 2 
eae 2 yx af | 
Ox, ~~ Age la Loe: | yy)? (1 m) x2 | yy” ” 
Pp Divan ; af 
= 14 SP leet (80) 
ve | ( me (x? cee i | . sa Abb. 5. Scheibe mit Einzellast. 
P. ZY" % x 
Auf den Randern y= +6 und y = — 6 wirkt demnach zunachst 
P ; 2.08 At (4 5 | | 
Ove Olas Oy=—6 — A\ (1 I m) (x? ae 2)? | ( 4 m) x? + 6? ? (30a) 
2 0? x x“ 
Ty =6 — Ty=—6 = ha Fla = m) (x? ae 62) (1 m) x ae =| : | 


1 K. Girkmann, Ing.-Arch. 11 (1940), Seddon 
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Mit (30a) wird dann, damit die Spannungen auf den Randern verschwinden, 


P DO ee 6 —ilx 
y Q3 = Ine | (1 +- m) (x? a 522 (3 ie m) x2 = | é A dx ’ 
ar +95= 0 Cy 
Pi TI 
U 2 07% x Sahih 
G2 + a= Ind | a Hm) Gap gap TE ™) aa sf ie 


Um die Integrale (31) auszuwerten, setzt man 


(orev be bale . lu ES. 
ae 0<7s 
leper =f (a2 2ye 


iG 


Die Integrale lassen sich mit x= u + iv itber die Residuen bal 
rechnen. Die Nullstellen der Nenner sind x)= +106, wobei x) 
in J, und J, ein Pol erster Ordnung, in J; und J, ein Pol zweiter; 


u 
= Ordnung ist. Mit dem durch Abb. 6 gegebenen Integrationsweg | 
fiir J < 0 ergibt sich | 
—|A|6 7 | 
= — — { — Sm [2| 6 
' J Lie, © Aad aI me 
Abb. 6. Integrationsweg " p 
zur Auswertung yon Ty4 a Ue <I [Al 5 Par TU i, £2 (Ns 
Js are x. ea 7 a O:|2|)\ e—1ahea 


Diese Ausdriicke setzt man in (31) ein und erhalt 


(a Nae Feta +m) ( m) (1 — 6|A|)e— 1418 — (3 +. m)e—l19] , | | 
(31a) 
(P2 +9) = ol (1 +m) Ade—l419-+-(1 — m)e—!419). 


Mit (31a) geht man in (22) und (27) ein und bekommt die Verschiebungen und Spannungen 


sty CO ae ee) 
— i thx u Cv} 6 Si 
re oe ee aes | (P1 t)eaeis a eae 
Beco cara) 
t tae iis Gindd Sindy , 
axe | mt elee tae: | +edaanecane Ora 
dk es ES 5 
eee. : the 1 Cof 25 CojA 
Vie ae (i ~ 5) Be Moe up| ow 1) Sate zagte* aa, 
= i x 1 : Si AOE A . | 
= fie + Pq) Lge'** dA tel ©: + Pa) Se ihe d), (32) 
+ | 
1 
Ox, = | 9) Kr + (Pp, +4) L,Ajei** dA, | 
S285 | 
} 
1 
aes (ee 2 bia — 3) 4? Ks + (P24 G4) A? Lele! ** di, 


vie AICI Pa aise aC eee 
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Auf dem Rand y= 6 greifen noch die Auflagekrifte an. Die Unterfunktion fiir eine Einzel- 
last (vgl. Abb. 7) ist 
a+é 
ae elon ls ey 1 - 
PA = — 7 lim po e—? do stoa tah evi At: (33) 


a—s 


Aus (33) erhalt man durch Grenziibergang fiir ein an der Stelle x= a 


angreifendes Moment Abb. 7. Zur Berechnung der 
Unterfunktion fiir eine 


iu=— = M,e—i*4, (34) Einzellast. 


Geht man mit (33) und (34) in (22) ein, so sind die Verschiebungen und Spannungen 


+ 2 +2 


peal z ws = : i — | CojAd Sindy Sin 12d Coj A 
asks ig i K Sok meal SojAd Sttry | SinArdGojry_ 
Ys rg | a Pi BalAe’ cds ao 
6) — oo x het4* dd), 
+ 20 + co 
1 = . 1 — | CojAdCoja SGinAd Sind 
eee A eS aK K idx 76 ey ys yi 
%3 aay roe ite) Aaa 2 & | eres ae Sly | 
— co — © x Aei** dd), 
i p (35) 
Ox3 = es 9, (K, + K,) A2et** dd, 
+ © 
1 = : 
ia a | eK + Kyteroan, 
+ 2 
t = : 
tT =— 55 | Get Kye eran, 
wobei 


= A : : L : : 
Caage (ensshe + et42) + + M, (e—i4a — eit a) 


ist. Damit ergeben sich die Spannungen 
Ce =x, 1 Cx, Cup Oy = Cy, Fy, + Oy B= YH 12 1-7 - 
Um die unbekannten Auflagekrafte und Einspannmomente zu ermitteln, geht man von 
ui =u, -+u, und v= v, +; aus. Die Auflagekrafte und Einspannmomente erbhalt man aus 


der Bedingung, daB die Integranden in u und v endliche Grenzwerte fiir 4 0 haben miissen. 
Es ergeben sich vier lineare Gleichungen zur Bestimmung der unbekannten Auflageki afte und 


Einspannmomente. 


(Eingegangen am 16. Dezember 1949.) 
Anschrift des Verfassers: Dr.-Ing. Hans Jung, Stuttgart-N, Technische Hochschule, Keplerstr. 10 
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Gasstrémung konstanten Querschnitts mit Warmezufuhr. 


Von O. Schrenk. 


1. Anschauliche Uberlegungen. Die Gasstrémung konstanten Querschnitts mit Warmezufuhr 


wurde neuerdings mehrfach theoretisch betrachtet!. Die dabei gefundenen Ergebnisse standen 
teilweise in einem scheinbaren oder wirklichen Widerspruch miteinander. 

Es ist indessen nicht schwer, das Problem auf anschaulich physikalischem Wege zu iiber- 
blicken, und auch einige Ergebnisse ohne Formeln abzuleiten. 

Die nachstehenden Uberlegungen zu dieser Frage* nehmen gewisse Analogien der Strémung 
in einer Heizkammer mit der in einer Lavaldiise zu Hilfe. 

Im Unterschallbereich der Lavaldiise nehmen Geschwindigkeit und Machsche Zahl zu, Druck, 
Temperatur und Dichte aber (adiabatisch) ab, wenn der Strémungsquerschnitt sich verengt. Die 
Energie der Masseneinheit bleibt dabei konstant. 

Will man nun in einer durchstrémten Kammer konstanten Querschnitts die Geschwindig- 
keit anwachsen lassen, so mu man hier das Gas durch Warmezufuhr ausdehnen. 

Die Temperatur T wachst dabei, im Gegensatz zum Fall der Lavaldiise, im allgemeinen in 
Strémungsrichtung an, abgesehen von einem noch zu erérternden Ausnahmegehiet. Der Druck p 
dagegen nimmt nach dem Impulssatz in Strémungsrichtung ebenfalls ab. 

Die Machsche Zahl 


M=~=—— = 
Co Veer, , 
(w Strémungsgeschwindigkeit, c Schallgeschwindigkeit, ~ Verhaltnis der spezif. Warmen c,/¢,, 
g Erdbeschleunigung, R Gaskonstante) wachst in der geheizten Parallelstrémung stets an, 
eleichviel ob dabei die Temperatur zu- oder abnimmt, da die Zunahme von w stets iiberwiegt. 
Diese an sich plausible Tatsache lat sich mit Hilfe der analytischen Beziehungen leicht be- 
statigen. 

Vom TemperatureinfluB abgesehen, besteht also eine Analogie zwischen der Lavaldiise und 
der zylindrischen Brennkammer; allerdings ist in dieser die Zustandsanderung nicht adiabatisch. 

Der auf den Temperaturverlauf zuriickgehende Unterschied besteht im wesentlichen darin, 
daB die Schallgeschwindigkeit mit der Temperatur steigt, so da gleiche Machsche Zahlen in der 
Heizkammer bei héheren Geschwindigkeiten erreicht werden als in der Lavaldiise, wenn man von 
gleichen Hintrittsbedingungen ausgeht. 

Die Analogie besteht auch hinsichtlich des besonderen Verhaltens der Strémung im Schall- 
querschnitt. 

Die Lavaldiise ist fiir M = 1 auf einer kurzen Wegstrecke zylindrisch wie die Heizkammer. 
Da nun die Lavaldiise keine Energiezufuhr besitzt, mu und darf auch der zylindrischen Heiz- 
kammer bei M= 1 keine Energie zugefithrt werden. 

Weiter sieht man entsprechend dem Bisherigen, daB eine Steigerung der Geschwindigkeit 
itber die Schallgeschwindigkeit hinaus Warmeabfuhr erfordert. Die Heizkammerstrémung be- 
sitzt also bei M= 1 ein Energiemaximum. 

Auch die erreichbare Héchsttemperatur besitzt einen Héchstwert. Dieser befindet sich aber 
bereits im Unterschallbereich, da bei kleinen Machschen Zahlen die Temperatur auf jeden Fall 
mit der zugefithrten Energie zunimmt und andererseits nach dem Vorausgegangenen bei M= 1 
ein negativer Temperaturgradient herrscht wie bei der Lavaldiise. 

Da8B bei M= 1 die zugefiihrte Energie ein Maximum hesitzt, scheint auf den ersten Blick 
itberraschend, da man ja auch bei Uberschallgeschwindigkeit weiteren Kraftstoff verbrennen kann. 
Und man kann auch, beispielsweise mit einem stromaufwarts befindlichen Kompressor von ge- 
niigender Leistung, eine beliebige groBe Férdermenge durch die Heizkammer pressen. 

Auch in dieser Frage schafft die Analogie mit einer in eine zylindrische Rohrleitung eingesetz- 
ten Diise Klarheit. Solange diese Diise nur einen konvergenten Teil besitzt, kann man in deren 
engstem Querschnitt ohne Riicksicht auf die Leistung des Kompressors nie die Schallgeschwin- 


: Schrifttum siehe J. V. Foa und G. Rudinger, J. Aer. Sci. 16 (1949), S. 84 und S. 566. 
Ecstmals vorgetragen im Juli 1949 in St. Louis im ElsaB. 
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digkeit itberschreiten. Die Wirkung einer Leistungssteigerung des Kompressors besteht nur 
darin, das die Luft vor dem Eintritt in die Diise komprimiert wird. Durch diese Kompression 
wird zwar die geférderte Masse erhéht, wie erwartet, aber das Niveau der Geschwindigkeiten und 
Machschen Zahlen vor und in der Diise in dem erforderlichen Mafe herabgesetzt. Man kann dies 
als den ,,Sperr- oder Drosseleffekt*: einer Diise bezeichnen. 

Der gleiche Sperreffekt tritt auch vor und in der Heizkammer auf, da auch hier nach dem 
Impulssatz, wie bereits gezeigt, ein Druckgefalle zwischen Ein- und Austritt herrscht. 

So sorgt der Strémungsmechanismus selbst automatisch dafiir, daB die Schallgeschwindigkeit 
in einer Heizkammer ohne Warmeabfuhr nicht iiberschritten wird. 


2. Rechnerische Zusammenhiinge. Bei der analytischen Durchrechnung komplizierterer Stré- 
-‘mungsanordnungen, die eine Heizkammer als Element enthalten, kommt man ohne besondere 
VorsichtsmaBnahmen ebenfalls ganz von selbst nur auf solche Liésungen, die mit dieser Grenze 
vereinbar sind, wenn man die im folgenden zusammengestellten Grundgleichungen (1) bis (4) 
oder eine entsprechende Anzahl von daraus abgeleiteten Gleichungen [z. B. (5) bis (9)] einfiihrt 
und die GréBen des Gas- und Strémungszustands an ihrem Eintritt und Austritt als Unbekannte 
behandelt. Zu diesen Gleichungen treten dann weitere hinzu, die den iibrigen Elementen der Ge- 
samtanordnung und den Randbedingungen des Problems entsprechen. 

Um Aussagen iiber die Brennkammer selbst zu machen, ist es aber nicht nétig diese tibrigen 
Gleichungen des Gesamtsystems heranzuziehen oder gar vom Strémungszustand ohne Heizung 
auszugehen!, 

Die nachstehende analytische Rechnung kann sich daher auf die eigentliche Brennkammer 
‘beschranken. 

Sie dient zur Bestadtigung und quantitativen Festlegung dessen, was im Vorausgegangenen 
schon auf Grund anschaulicher Uberlegungen gefunden worden ist. 

Der Zeiger 1 bedeutet im folgenden den Eintrittsquerschnitt der Brennkammer; die G1réBen 
ohne Zeiger gehéren zu einem beliebigen Schnitt innerhalb oder am Austritt der Kammer. Die 
Wandreibung ist vernachlassigt. Die Energiezufuhr ist stetig und in strémungssenkrechten 
Schnitten konstant angenommen. 

Man hat dann, wenn noch g die Gasdichte bezeichnet, die bekannten Grundgleichungen 


ow = 0,w, (Kontinuitat), (1) 
P — Pi_ (Gaseleichung), 2 
ey re et ( : 8) (2) 
p +ew = p, +e; (Impulssalz), (3) 
q_ 8pm , wi 
In der letzten Gleichung ist q die der Masseneinheit vom Kintrittsquerschnitt ab zugefiihrte 
Warmeenergie in a , A das mechanische Warmedquivalent, c,,, die im voraus leicht 
g se 


abzuschatzende spezifische Warme zwischen T und T,, die etwas gréBer als cp, ist. 
Wenn man das Aufheizungsverhaltnis 


und das Entspannungsverhaltnis 


einfiihrt, erhalt man aus (1) bis (3) eine Reihe von Beziehungen, die man als Funktionen der vier 


Verdnderlichen 
t, yp, M, und M 


ausdriicken kann, wobei jeweils zwei dieser Veranderlichen die beiden anderen bestimmen. 
Zuniachst findet man aus (3) allein einen anderen, oft niitzlichen Ausdruck fiir den Impulssatz, 


wenn man 
xP =e 
Q 
1 Foa und Rudinger, die diesen Weg gehen, sind iibrigens in einigen Anwendungsbeispielen zu zahlen- 
mafBig identischen Ergebnissen gelangt wie der Verfasser. 
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einfiihrt, namlich 


1 +x, M3 
VS ileeioeViaee 
Weiter ergibt sich aus (1) bis (3) 
M, 
O4 
y=140 
t=2 
g3P .| 
3 
4 
G2 & 


01\— 


a2 40 


zeichneten Werte: 1 a 
M = — (= 0,845 fir % = 

Vz 
(M ist unabhangig vom Zustand im Lintrittsquerschnitt), 


ae 


* Unter Beschrinkung auf den rechten Ast der Kurven, s 


x» M? (1+ ~,M?2)? 


~ 4,3 (1 - x M2? ’ (6) | 
; es 
mMi—S=OP, (| 
(Vp) x i 

Me ee 8) | 

i («—y)¥ (8) | 


Ferner gewinnt man aus (1) und 


(2) unmittelbar die einfache Be- | 


ziehung 
Wixcosh Fe 


ee Gee (9) | 


Wy yp 


Die Gleichungen (5) bis (9) lassen | 


sich fiir die Anwendung noch in ver- 


schiedener Weise umformen und | 


kombinieren. 


In Abb. 1 ist M, nach (7) far | 
x, = 1,40 abhangig von py und ¢ dar- | 
gestellt. Man sieht hier die ,,Sperr- | 


wirkung‘* der geheizten Kammer: 
pw ist bei gleichem M,, und ebenso 


M, bei gleichem y, um so kleiner je | 


gréBer t ist. Auch der bereits physi- 
kalisch erérterte Héichstwert der 


Aufheizung fiir ein gegebenes M, er- | 


scheint in dieser Abbildung (Scheitel- 
punkte der Kurven), und zwar ist 
er um so kleiner, je gréBer M, ist. 

DaB diese Héchsttemperatur bei 
einer lokalen Machschen Zahl M< 1 
liegt, bestatigt die Abb. 2, in der M 
abhangig von y und ft dargestellt ist. 

Diese Abbildung zeigt auBerdem, 
daB die kleinen Werte von yp nur er- 
reicht werden, wenn man M > 1 zu- 
laBt, das heiBt, wie schon physi- 


kalisch gezeigt wurde, mit Warme- | 


abfuhr. Ob Warme zu- oder ab- | 


gefiihrt werden mu8, dariiber kann 


allerdings diese Abbildung ebenso- — 


wenig aussagen wie die Gleichungen 
(5) bis (9), da die Energiegleichung 
dabei nicht benutzt wurde. 

Fiir das unterhalb M= 1 liegende 


, temperaturmaximum gibt die Rech- 
nung, wie bekannt, die nachstehen- | 
den, durch einen Querstrich gekenn- | 


1,40) 


. die folgenden Ausfiihrungen. 
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p= (Vz —-yz—1) 


und 
s2GaaMm 
4 2, M? 
Fir die Schallgrenze (Zeiger *) mit 
i he | 
erhalt man in entsprechender Weise 
1+ x, M2 
eee at 1 
geen TU Zoe 
ne 2 
5 = pgs (= 0.834 fir x* = 1,40), 
== aes yrlae 0,973 fiir x* = 1,40). 


Die bisher nicht benutzte Energiegleichung (4) gibt Aussagen iiber das Energiemaximum und 
die Verteilung der zugefiihrten Warmeenergie auf den rein thermodynamischen Anteil (Warme- 
inhalt) und die kinetische Energie. Man kann 


(4) in der dimensionslosen Form i ae 

A i Z Bom a 
SS Vad Pomme 7 7a alee 
& opm Ty Caged hal Cpm : a ( 1 


schreiben; diese Funktion ist in Abb.3 fiir 
den Sonderfall M,= 0,2 als Beispiel dar- 
gestellt. Man findet hier die zugefiihrte Ge- 
samtenergie q/gc,,, 7, und den rein thermo- 
dynamischen Anteil r— 1 abhangig von yp 
aufgetragen!. Der schraffierte Teil zwischen 
den beiden Kurven bedeutet die kinetische 
Energie. AuSerdem ist M als Funktion von 
py eingetragen. Man sieht, daf jenseits y und 


M die in kinetische Energie umgewandelte 
Energie gréfer ist als die zugefiihrte, so daB 
die Gastemperatur hier trotz der Heizung 
abnimmt. 

Die Darstellung der Abb. 3 1aBt sich ohne 
Schwierigkeit verallgemeinern und von M, 
unabhangig machen, doch ergeben sich dabei Aue. 
keine wesentlichen neuen Gesichtspunkte. 


3. Beziehungen zum Verdichtungssto8. Man kann nunmehr auch eine gewisse Analogie 
zwischen der Strémung mit Warmezu- und -abfuhr und der des geraden Verdichtungsstofes her- 
stellen. Die drei Gleichungen (1) bis (3) gelten in genau gleicher Form auch fiir den geraden 
StoB. Bei diesem tritt nur an die Stelle von (4) die Bedingung qg= 0. Folglich definieren jeweils 
zwei Punkte gleicher Energie der Abb. 3, einer rechts und einer links vom Energiemaximum 
(M* = 1), einen Verdichtungssto8. Wahrend aber im Falle des StoBes (q= 0) aus Entropie- 
griinden nur die unstetige Verdichtung vom Uberschall- auf den Unterschallzustand méglich ist, 
kann mit Hilfe der unbeschrankten Warme- und Kaltequellen der jetzigen Anordnung eine 
stetige Zustandsdnderung stattfinden. . 

Nach den bisherigen Betrachtungen kann man sich schliefSlich auch ohne Schwierigkeit quali- 
tativ dariiber Rechenschaft geben, was geschieht, wenn die Gas- und Strémungszustande im 
umgekehrten Sinne durchlaufen werden. 


1 Um die iibliche Vorstelluug einer Strémung von links nach rechts zu behalten, ist der Anfangs- 
punkt y=1, d.h. der Eintrittsquerschnitt ans linke Ende der Abszissenachse verlegt. Die Abbildung 
umfaBt die Darstellung aller Heizkammern mit M,=0,2; bei y=1 liegt stets der Eintritt, der Aus- 
tritt kann bei beliebigem y gewahlt werden. Der y-MaQstab ist gegeniiber dem Langenmafistab der 
Kammer natiirlich in einer nicht genauer bekannten Weise verzerrt. 
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Im reinen Unterschallbereich bietet sich dabei kein neues Problem: durch Warmeentzug; 
wird das Gas verdichtet und die Strémung verlangsamt. Theoretisch besteht diese Méglichkeitt 
bis zur Erreichung der Geschwindigkeit null bei T= 0. Der Druck erreicht in diesem Grenzfall| 
einen endlichen Maximalwert gleich dem Doppelten des Wertes beim Temperaturmaximum.. 

Fiihrt man einer Uberschallstrémung Warme zu, so muf man zwei Falle unterscheiden: 

Wenn die Warmezufuhr nicht das durch das Energiemaximum (Abb. 3) gegebene Maf iiber-- 
schreitet, so wird die Strémung verzégert, bleibt aber iiberkritisch. Eine Riickwirkung der Hei-- 
zung auf stromaufwarts liegende Bereiche, also eine ,,Sperrwirkung“‘ wie im Unterschallfall, ist 
nicht vorhanden. 

Fiihrt man aber nach Erreichung der Schallgeschwindigkeit noch mehr Energie zu, dann ist 
die bisher vorausgesetzte und im Beispiel der Abb. 3 dargestellte Strémungsform nicht mehr 
existenzfahig. Es wird jetzt ein VerdichtungsstoB auftreten, der im Schallquerschnitt mit der; 
Intensitat null beginnt, um bei weiterer Energiezufuhr stromaufwarts in die Heizkammer hinein: | 
zuwandern. Man weif ja, daB eine Verdichtungswelle endlicher Intensitat gegen die Uberschall-: 
strémung stromaufwarts vordringen kann, da ihre Wellengeschwindigkeit gréBer als die Schall- 
geschwindigkeit ist. Die Umwandlung eines Teiles der Uberschallstrémung in eine Unterschall-- 
stromung erméglicht die Unterbringung gréBerer zugefiihrter Energiemengen. Der StoB wird 
sich dabei so einstellen, daB der Enddruck der Brennkammer sich der Austritts-Druckbedingung; 
anpaBt. 

Auf der Grundlage dieser, nur skizzierten, physikalischen Vorstellung kann man im Bedarfs-- 
falle ein quantitatives theoretisches Rechenverfahren entwickeln. Doch ist diese Aufgabe, wie> 
es scheint, von geringerer technischer Bedeutung. 


(Eingegangen am 12. Januar 1950.) 
Anschrift des Verfassers: Prof. Dr.-Ing. 0. Schrenk, Paris XV, 403 Rue de Vaugirard. 
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Das in Wandnihe giiltige Geschwindigkeitsgesetz 
turbulenter Strémungen. 


Von J. Rotta. 


Bei turbulenten Strémungen an freien Wanden und in Rohren oder Kanalen hat sich gezeigt, 
da8 in Wandnihe ein universelles Geschwindigkeitsgesetz gilt, das die Form 


eles aie 
a ( . { (1) 
hat, wabei U die mittlere Strémungsgeschwindigkeit im Abstand y von der Wand, 


v*¥ = r/o (2) 
die mit der Wandschubspannung ft, und der Dichte @ gebildete Schubspannungsgeschwindigkeit 
und » die kinematische Zahigkeit sind. Dieser universelle Zusammenhang wurde zuerst von 
Prandtl fiir Rohrstémungen und Grenzschichten an Wanden bei konstantem AuBendruck nach- 
gewiesen, wobei die Wande hinreichend glatt waren. Neuerdings ist H. Ludwieg und W. Till- 
mann* der experimentelle Nachweis dieses Gesetzes auch fiir Grenzschichten bei beliebigem 
Druckgefalle gelungen. 


Das Gesetz (1) ist nun aber von der geometrischen Beschaffenheit der Wand, d. h. also von 
der Wandrauhigkeit abhangig. Insofern ist die Bezeichnung ,,universelles. Geschwindigkeits- 
gesetz‘‘ nicht ganz treffend; sie soll lediglich zum Ausdruck bringen, da das Geschwindigkeits- 
gesetz unabhangig von den duferen Strémungsbedingungen in groBem Wandabstand ist. Ein 
~universelles Gesetz in diesem Sinn kann natiirlich nur in einer sehr diinnen, wandnahen Schicht 
bestehen, die einen kleinen Bruchteil der gesamten Grenzschichtdicke bzw. — bei Rohrstré- 
mungen — der Rohrabmessungen ausmacht. Was dem Gesetz (1) aber die groBe Bedeutung gibt, 
ist die Tatsache, daB sich innerhalb dieser diinnen Schicht die gesamte Wechselwirkung zwischen 
der von der Zahigkeit erzeugten Schubspannung und der durch die turbulenten Geschwindig- 
keitsschwankungen bewirkten Schubspannung abspielt. Am auferen Rande dieser Schicht hat 
die Turbulenz bereits eine solche Intensitat, daB die Zahigkeitsspannung daneben vernachlassig- 
bar ist. Hieraus ergeben sich z. B. wichtige SchluBfolgerungen fiir die Behandlung turbulenter 
Grenzschichten, wie in einer anderen Arbeit von mir? gezeigt wird. 

Fir die Darstellung des Gesetzes (1) hat sich die Benutzung des Mischungswegansatzes von 
L. Prandtl* als zweckmaBig erwiesen. Die sich hiermit fiir glatte und rauhe Wande ergebenden 
Zusammenhange sind bekannt. Es laBt sich nun zeigen, daB mit dem Mischungswegansatz bei 
physikalisch sinnvoll gewahlten Randbedingungen auch der Ubergang in die laminare Unter- 
schicht verniinftig erfaBt wird und die Wirkung der Wandrauhigkeit anschaulich erklart werden 
kann. 

Nimmt man in dem Prandtlschen Ansatz fiir die Schubspannung 


iG qu | se (3) 


dy dy 
eine Linearbeziehung zwischen dem Mischungsweg | und dem Wandabstand y an und setzt fer- 
ner die Schubspannung 7 konstant gleich der Wandschubspannung 1), was unter den genannten 
-Voraussetzungen fiir die Giiltigkeit des universellen Gesetzes berechtigt ist, so ergibt die Inte- 
gration der Differentialgleichung (3) eine Naherung fiir den Geschwindigkeitsverlauf, die durch 
Messungen als brauchbar bestatigt wird. Innerhalb der laminaren Unterschicht, deren Dicke 
Oz sei, ist keine Turbulenz vorhanden, deshalb ist in (3) /= 0 zu setzen. Mit U= 0 fir y= 0 


1 L. Prandtl, Zur turbulenten Strémung in Robren und langs Platten, Ergebnisse der Aerodynamischen 
Versuchsanstalt zu Géttingen, S. 18. IV. Lfg. Oldenbourg. Miinchen 1932. 

2 H. Ludwig und W. Tillmann, Ing.-Arch. 17 (1949), S. 288. 

3 J. Rotta, Uber die Theorie der turbulenten Grenzschichten. Mitteilungen aus dem Max-Planck- 
Institut f. Strémungsforschung Nr. 1, Géttingen 1950. ' ‘ 

4 L. Prandtl, Fihrer durch die Strémungslehre, 3. Aufl. Braunschweig 1949. 
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folgt somit fiir 


*2 
VOL U= y=" 9 (Aa) 
Fiir den anschlieBenden, turbulenten Teil kann man ‘ 
= #(y — 02) (5) 


setzen!, wobei x als eine universelle Konstante angesehen wird, die etwa die GréBe x ~ 0,4 hat. | 
Der Anschlu8 an (4a) fiir y= 6; wird durch die Bedingung 


v6 
Use === m 


gesichert. Die Integration von (3) liefert ,wenn man / fiir y mit Hilfe von (5) substituiert, fiir 


pahcumelseft EPH] tm! 2/ ET] +4. | 


@ P| 

‘wh Abb. 1 zeigt (4a) und (4b) fiir v*0,/y = 6,7 | 

Pla und x—0,4im Vergleich mit Messungen von | 
H. Reichardt? an als hydrodynamisch glatt 
20 
yg so 
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Abb. 1. Geschwindigkeitsprofil nach (4a) und (4b) fiir die Abb. 2. Einflu8 schwacher Wandrauhigkeit auf die Dicke dy 
glatte Wand mit v* 07 /y = 6,7, x = 0,4, Versuchspunkte der laminaren Unterschicht nach Versuchsergebnissen 
nach H, Reichardt. von J. Nikuradse. 


anzusehenden Wanden. Fiir groBe y-Werte hat (4b) die bekannte asymptotische Form 


sella las ce \ 

mit 
1 v* dp 

C—— (In 4% — 1) apisre (7) | 


Bei rauhen Wanden sind die Verhaltnisse beziiglich der Randbedingungen weniger iibersicht- | 
lich. Man wird auch hier eine Bezugsebene fiir den Wandabstand y definieren, in der die Ge-: 
schwindigkeit U im Mittel verschwindet. Die Lage dieser Ebene, die sich vermutlich nur wenig | 
unterhalb der gré8ten Erhebungen befindet, ist nicht genau bekannt; sie ist aber hier auch nicht 
so sehr wichtig. Eine kleine Wandrauhigkeit hat nun zur Folge, daB die Strémung in der lami- } 
naren Schicht wellig wird; dieses begiinstigt die Entstehung der Turbulenz, Praktisch wirkt sich) 
dies in einer Abnahme der Dicke 6, der laminaren Unterschicht aus. In Abb. 2 ist die diesheziig- 
liche Auswertung der MeBergebnisse von J. Nikuradse? an Rohren mit sandrauhen Wanden\ 
wiedergegeben; k ist die Korngré®e der Rauhigkeiten. Die Dicke dz fallt mit wachsendem k ab 
bis sich etwa fiir v*k/y = 54 der Wert 6; = 0 ergibt. 


2) 


* Diese Annahme wurde auch in einer Arbeit von H. B. Squire, Phil. Mag. (7) 39 (1948), S. 7, gemacht;} 
jedoch wurde dort ein von (3) abweichender Ausdruck fiir die Schubspannung benutzt. Der berechnete | 
Geschwindigkeitsverlauf weicht nur wenig von dem hier erhaltenen ab. ae | 

* H. Reichardt, Z. angew. Math. Mech. 20 (1940), S. 297. 

* J. Nikuradse, Strémungsgesetze in rauhen Rohren. VDI-Forschungsheft 361 (1933). 
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Fiir noch gréRere Rauhigkeiten besteht dann keine laminare Unterschicht mehr, sondern die 
Turbulenzintensitat hat schon in der Ebene y = 0, in welcher der Mittelwert der Geschwindigkeit 


U verschwindet, eine endliche Gri- 
Be. Bei sehr groBen Rauhigkeiten 
ist diese Intensitat bei y= 0 schlieB- 
lich so stark, daf die Zahigkeit iiber- 
haupt keinen Einflu8 mehr hat. Man 


18\—— 


16 


Do 
vlog 


darf wohl] annehmen, da8 in diesem 
Fall an den einzelnen Bergen und 


Talern der Wand unmittelbar Ab- 


lésung der Strémung eintritt. In (3) “ 
wird dieser Sachverhalt dadurch aus- 
gedriickt, daB® der Mischungsweg fiir 


y = 0 eine endliche, von der Rauhig- 7% 
keit abhangige GréBe |, hat. An 
Stelle von (5) ist dann 


=htxy (8) ” 
zu setzen. Hiermit ergibt die Inte- 
gration von (3) mit U—0 fir y= 0 


van el vere), 


* 


ancl pany el 
xv*l, | 2 V y . g 


| 
jk ~ +V(S) 4 | 
x v* I, / a wee 
(eg \( v "A 
Nach den zitierten Messungen von J. 9 
Nikuradse an sandrauhen Rohren 
ist, wie Abb. 3 zeigt, der Anfangs- 
wert J, von k linear abhangig; fiir v*k/v 
=> 54 gilt 


OG v* i 


0,76 4 


(10) 


Auch die Gleichung (9) geht fiir groBe 
y-Werte in die Form (6) itber, wobei jetzt 


v i 
Hv*ly| 2 


eee 
ee eS 


ist. Fihrt man hierin v*l,/y nach (10) 
ein, so wird man fiir sehr groke v*k/y- 
Werte auf die von L. Prandtl angegebene 
Beziehung 


u=v(Lind +c) 


i 
C= —(In4x—1) 


(11) 


(12) 


gefithrt, wobei sich in unserem Fall C,— 8,37 ergibt. 


Zahigkeit. 
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Abb. 3. Anfangswert fiir den Mischungsweg bei starker Wandrauhigkeit 
nach Versuchsergebnissen von J. Nikuradse. 
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Abb. 4. Berechnete Geschwindigkeitsprofile fiir verschicden 


rauhe Wiainde mit x = 0,4. 


Die Form ist ganz unabhangig von der 
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Zum AbschluB dieser Ausfiihrungen sind in Abb. 4 die aus (4a), (4b) und (9) errechneten Ge- 
schwindigkeitsgesetze fiir verschieden rauhe Wande zum Vergleich miteinander aufgetragen. 
Das vorliegende Versuchsmaterial iiber Strémungen an rauhen Wanden reicht leider nicht aus, — 
um den dargestellten Geschwindigkeitsverlauf bis in unmittelbare Wandnahe mit Messungen ver- 
gleichen und die Lage der Bezugsebene y = 0 bestimmen zu kénnen. 


(Eingegangen am 12; Januar 1950.) 


Anschrift des Verfassers: J. Rotta, Gottingen (Hannovy.), Béttingerstr. 6/8. 
Max-Planck-Institut fiir Strémungsforschung Gottingen. 
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Vor kurzem erschienen: 


- Eimfuhrung in die Akustik 
Von 
Dr. phil. Ferdinand Trendelenburg 


Ve Honorar-Professor an der Universitat Freiburg i. Br. 


Zweite, umgearheitet Auflage } 
Mit 280 Abbildungen, VIII, 378 Seiten. 1950. Ganzleinen DMark 39.— 


Inhaltstibersicht: { 


1. Grundiegende Fragen der Schwingungstehre und der Wellenlehre: Einfache Schwingungen. Zusammengesetzte 
Schwingungen. Fourierdarstellung, Lissajousschwingungen, Freie und erzwungene Schwingungen (Systeme. von 
einem Freiheitsgrad). Schwingungen von Systemen mit nichtlinearen Eigenschaften. Koppelungsschwingungen 
(Systeme von mehreren Freiheitsgraden). Selbsterregte Schwingungen. Wellengleichuny. Die verschiedenen Wellen- 
arten. Eigenschwingungen yon Luftsdulen, Saiten, Membranen, Staben. Hinflu8 einer Bewegung von Schallquelle 
oder Schallempfanger. Dopplereffekt. — I]. Schallfeldgré6en und ihre Messung: Schalldruck, Bewegung, Schnelle, 
Temperaturschwankung, Dichteschwankung. Frequenz, Wellenlinge, Schallgeschwindigkeit. Schallstarke, Schall- 
leistung. Lautstarke. — [Il./Schallerzeugung: Grundlegende theoretische Bemerkungen zur Schallabstrahlung. 
Mechanisehe Schallsender, Musikinstrumente. Die menschliche Stimme, Elektrische Schallsender. Thermische 
Schallsender. — IV. Schallausbreitung: Schallgeschwindigkeit. Huygens’sches Prinzip, Reflexion, Beugung, Brechung, 
Schallabsorption. Vorgange, in geschlossenen akustischen Systemen (Resonatoren, Filter, Leitungen). Raum-\und 
Bauakustik. — WV. Schallempfang und Schallaufzeichnung: Wirkungsweise und Bauart technischer Schallempfanger. 
Gerichteter Schallempfang. Eichung von Schallempfanger. Das Ohr als Schallempfanger. Schallaufzeichnung. — 
Vi. Schallanalyse. Physikalische Eigenschaften natiirlicher Schallvorgange: Verfahren zur Schallanalyse. Physka- 


lische Eigenschaften natiirlicher Schallvorginge. — VII. Anhang: Benennungen in der Akustik. Zusammenstellung 


praktisch wichtiger akustischer Formeln. — Sachverzeichnis. » 


Fir die zweite Auflage der >, Einfihrung in die Akustik** wurde die Einteilung des Stoffes beibehalten. Am ersten» 
die grundlegenden Fragen der Schwingungs- und Wellenlehre behandelnden Abschnitt, brauchte naturgem48 nur 
wenig geindert zu werden. Alle ibrigen Abschnitte des Buches wurden im Hinblick auf die schnelle Weiterentwick- 
lung, die die Akustik seit Erscheinen der ersten Auflage genommen hat, weitgehend umgearbeitet und erginzt. Eine 
Reihe von Abschnitten wurde praktisch vollig neu geschrieben. | 


' 


. 


Der Kreisel 
Seine Theorie und seine Anwendungen 
Von 
Dr. R, Grammel 


o. Professor an der Technischen Hochschule Stuttgart 


Zweite, neubearbeitete Auflage 
Erster Band: 
Die Theorie des Kreisels 


Mit 137 Abbildungen. XI, 281 Seiten. 1950. DMark 30,—; Ganzleinen DMark 33.— 


j Inhaltstibersicht: 
Finleitung. — Grundlagen: Grundlagen derVektorrechnung. ~~ Grundlagen der Mechanik. — Der Tragheitstensor, — 
Der symmetrische Kreisel: Der kraftefreie symmetrische Kreisel. — Die gefihrte Bewegung des symmetrischen 
Kreisels. — Der symmetrische Kreisel unter Zwang und Sto8. — Der schwere symmetrische Kreisel. — Der Einfluf 
der Reibung. — Der unsymmetrisshe Kreisel: Der kraftefreie unsymmetrische Kreisel.—-- Der schwere unsymmetrische 
Kreisel. — Besondere Probleme: Kreisel im erweiterten Sinne. Gyroskopische Systeme. — Darstellung der Kreisel- 
bewegungen durch Thetafunktionen. — Namen- und Sachverzeichnis. : 


: 


Im September erscheint: } 
Zweiter Band: 


Die Anwendungen des Kzeiseis 
Mit 133 Abbildungen. Etwa 240 Seiten. DMark 30.—; Ganzleinen DMark 33.— 


Inkhaltsitbersicht: 


Kreiselwirkung bej Radsatzen. Kollermihlen. KritischeDrehzahlen von Rotoren. Fahrzeuge (Bahnen, Kraftfahrzeuge, 
Schiffe, Zweirad). Flugzeuge. — Kreiselgerate. KompaBkreisel (Gyroskop, Inklinations- und Deklinationskreisel). 
Kreiselkompasse (EinkreiselkompaB, Mehrkreiselkompasse). Kiinstliche Horizonte (Kreiselpendel, Pendelkreisel). 
Wende- und Lagekreisel (Wendezeiger, Kurskreisel, Richtkreisel, Stitzkreisel). Sonstige Kreiselgerate (Differentiir- 
und Integrierkreisel, Reglerkreisel). — Unmittelbare Stabilisatoren. Richtkreise! (Erde, geworfene Kérper), Stitz- 
kreisel (Einschienenbahnen, Geradlaufer). Diampfkreisel (Schiffskreisel). : 
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Vor kurzem erschienen: 


Integraltafeln 
Sammlung unbestimmter Integrale elementarer Funktionen 
Von 


Dr.-Ing. W. Meyer zur Capellen 
Aachen ; 


VIII, 292 Seiten. 1950. Ganzleinen DMark 36.— 


Inhaltstibersicht: 

Vorbemerkungen: Regeln zum Integrieren. — Regeln zum Differenzieren. — Weitere Hilfsmittel zur Integration, — 
Zum Gebrauch der Tafeln. — Integrale algebraischer Funktionen: Integrale rationaler Funktionen. — Integrale 
irrationaler Funktionen. Integrale algebraischer Funktionen, die auf elliptische Integrale fahren. — Integraletranszen- 
denter Funktionen: Exponentialfunktion und Logarithmus. — Integrale trigonometrischer und zyklometrischer Funk- 
tionen. — Hyperbel- und req-Funktionen. — Produkte algebraischer und transzendenter Funktionen: Exponential- 
funktion und Logarithmus. — Trigonometrische und zyklometrische Funktionen, — Hyperbel- und Wrea-Funktionen. 
— Produkte transzendenter Funktionen untereinander: Integrale von der Form fg(x) In xd«. —- Integrale von der 
Form fexg(x)dx. —’Zusammenstellung einiger wichtiger Konstanten, Reihen und Funktionen. — Schrifttum. 
Fiir die Auswahl der gebrachten Integrale — es diirften rund 3000 sein —- war maBgebend, daB der Benutzer einer- 
. seits fertige Integrale findet, andererseits auch Rekursions- oder Hilfsformeln, mit denen er weiterarbeiten kann. 


Solche Formeln wurden aber trotzdem fiir eingelne Sonderfalle ausgewertet, damit der Benutzer die Anwendung 
kennenlernt. 


Einteilung und Aufban erfolgten ausschlieBlich aus dem Gesichtspunkt des praktischen Gebrauches, d. h. der Be- 


nutzer, auch wenn er mathematisch nicht sehr geschult ist, soll méglichst schnell das gesuchte Integralfinden. Daher 
wurde, wie die Uberschriften und das Inhaltsverzeichnis zeigen, sehr weitgehend unterteilt, und zwar so, daB die 
Uberschrift einen Hinweis auf den Ort des Integrals gibt. — Zum Differenzieren. kann die Tafel ebenfalls benutzt 
werden, die Formeln sind dann nur von rechts nach links zu lesen. 


Das kurze Schrifttumsyerzeichnis enthalt die wichtigste benutzte und die im Hauptteil angezogene Literatur. 


Matrizen 


Eine Darstellung fiir Ingenieure 
‘Von 
Dr.-Ing. Rudolf Zurmihl 


Mit 25 Abbildungen. XV, 427 Seiten. 1950. Ganzleinen DMark 25.50 


Inhaltsibersicht: 


1, Einleitung. — I. Kapitel: Der Matrizenkalkiil: 2. Grundbegriffe und einfache Rechenregeln. 3. Matrizen und 
Vektoren. 4. Matrizenmultiplikation. 5. Kehrmatrix und Matrizendivision. 6. Lineare Transformationen. 7. Ortho- 
gonale Transformation. — II. Kapitel: Der Rang: 8. Determinanten. 9, Lineare Abhingigkeit und Rang. 10. Theorie 
der linearen Gleichungen. 11. Aquivalenz und Rangbestimmung. — III. Kapitel: Formen und Transformationen: 
12. Bilineare und quadratische Formen, 13. Koordinatentransformationen, — IV. Kapitel: Das Eigenwertproblem: 
14. Charakteristische Zahlen und Eigenvektoren, 15. Symmetrische Matrizen. 16. Allgemeinere Eigenwertprobleme. 
17. Komplexe Matrizen. — V. Kapitel: Struktur der Matrix: 18. Die Minimumgleichung: 19, Elementarteiler, Klassi- 
fikation. 20. Die Normalform, 21. Hauptvektoren. Transformation auf Normalform. 22. Matrizenfunktionen und 
Matrizengleichungen. — VI. Kapitel: Numerische Verfahren: 23. Auflésung linearer Gleichungssystemedurch Matrizen- 
multiplikation, 24. Iterative Behandlung linearer Gleichungssysteme. 25. Iterative Bestimmung der gréBten charak- 
teristischen Zahl. 26. Bestimmung héherer Eigenwerte. VII. Kapitel: Anwendungen: 27. Matrizen in der Elektro- 
technik. 28. Matrizen in der Schwingungstechnik. 29, Systeme linearer Differentialgleichungen. 30. Differential- 
matrizen und nichtlineareTransformatione n.31.Tensoren.32. Matrizen in der Ausgleichsrechnung, — Sachverzeichnis. 


Die Lehre von den Matrizen ist eine Lehre von den linearen Beziehungen. Diese aber spielen in der gesamten mathe- 
matischen Naturbeschreibung und dementsprechend auch in der Ingenieurwissenschaft eine hervorragende’ Rolle. 
Es gibt kaum ein der Rechnung zugingliches Gebiet in Physik und Technik, in dem man nicht mehr oder weniger 
zwangsliufig auf Bezichungen linearer Art gefiihrt wird. 


In den Matrizen hat sich die lineare Algebra ein Ausdrucksmittel von ganz besonderer Anpassungsfahigkeit geschaffen. 
Mit ihrer Hilfe werden einheitlich faBbare, in der. gewéhnlichen Formelsprache jedoch nur schwerfiallig darstellbare 
Operationen durch Formeln yon uniibertrefflicher Kiirze und Sinnfalligkeit wiedergegeben, die den Blick stets auf 
das Wesentliche lenken. Die Matrizenrechnung ist daher in besonderer Weise geeignet, die so vielgestaltige Welt 
der linearen Transformationen auch dem Nichtmathematiker und Ingenieur zu erschlieGen. 


Das vorliegende Buch wendet sich demgemifs nach Darstellungsart, Aufbau und Stoffauswahl an den Ingenieur, 
und zwar ganz besonders auch an jenen, der nicht tiber eine mathematische Sonderausbildung yerfiigt.' * 
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